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Abstract:  Complaints  about  Army  training  noise,  especially  artillery  noise,  are 
increasing  and  are  impacting  soldier  training  opportunities.  One  suggested  mitiga¬ 
tion  method  is  to  use  a  gravel  pad  near  the  noise  source  to  reduce  blast  noise. 
Measurements  were  conducted  to  assess  this  method  by  detonating  C4  charges 
located  over  a  15-  x  15-  x  1.5-m-thick  gravel  pad  or  over  undisturbed  ground  and 
recording  the  acoustic  and  seismic  waveforms  at  various  distances  from  the 
source.  The  measurements  recorded  by  ERDC-CRREL  personnel  at  propagation 
distances  from  10  to  400  m  are  documented  in  this  report.  Additional  reports 
documenting  the  longer  distance  measurements  and  analyzing  the  measurements 
are  planned. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes.  Citation  of 
trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products.  All  product  names  and 
trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to  be  construed  as  an  official  De¬ 
partment  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 

DESTROY  THIS  REPORT  WHEN  NO  LONGER  NEEDED.  DO  NOT  RETURN  IT  TO  THE  ORIGINATOR. 
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Preface 

This  report  was  prepared  by  Frank  E.  Perron  Jr.,  Stephen  N.  Decato,  Dr.  Donald 
G.  Albert,  and  David  L.  Carbee,  Geophysical  Sciences  Branch,  Cold  Regions  Re¬ 
search  and  Engineering  Laboratory,  U.S.  Army  Engineer  Research  and  Develop¬ 
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Many  personnel  contributed  to  the  success  of  this  test  program.  Gary  Leadore 
(ATC)  was  the  Test  Director,  Danny  Deal,  Range  Manager,  and  Martha  Tum- 
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Samantha  Rawlings.  ERDC-CRREL  participants  were  Don  Albert,  Frank  Perron, 
Steve  Decato,  and  Dave  Carbee.  The  Netherlands  participants  were  Frank  van  den 
Berg,  Jaap  van’t  Hof,  and  Frits  Van  der  Eerden,  TNO,  and  Erik  van  Arkel,  MOD- 
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Sgt.  Kone,  and  SPC  Hamlin,  site  support  personnel.  The  authors  thank  all  of  them 
for  assistance  during  the  tests,  and  they  thank  Dr.  Larry  Pater  of  CERL  for  fund¬ 
ing  their  participation  in  this  study. 
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Wuebben,  Acting  Director,  CRREL.  The  Commander  and  Executive  Director  of 
the  Engineer  Research  and  Development  Center  is  COL  Richard  B.  Jenkins.  The 
Director  is  Dr.  James  R.  Houston. 
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1  Introduction 

The  growth  of  civilian  populations  near  Army  training  facilities  has  led  to  increas¬ 
ing  complaints  about  Army  training  noise,  especially  the  noise  produced  by  artil¬ 
lery  fire.  Many  approaches  are  being  considered  to  reduce  the  noise  complaints. 
This  report  documents  the  seismic  and  acoustic  measurements  made  by  CRREL 
personnel  during  a  joint  test  of  one  possible  approach,  the  use  of  a  gravel  pad  to 
reduce  blast  noise.  These  measurements  were  conducted  to  determine  the  absorb¬ 
ing  effect  of  a  15-  by  15-  by  1.5-m-thick  gravel  pad  (Figure  1)  on  detonated  single 
brick  Ml  12  C-4  charges.  The  tests  were  conducted  on  the  TW  Range  at  the  Aber¬ 
deen  Test  Center  (ATC),  located  in  Aberdeen,  MD,  from  13  to  17  June  2005. 
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2  Approach 

This  test  was  to  determine  whether  the  blast  absorber  approach  can  provide  a  sig¬ 
nificant  reduction  (at  least  3  dB)  in  blast  noise  at  realistic  distances  from  the 
source.  If  successful,  further  development  of  this  method  is  planned.  The  concept 
of  the  blast  absorber  approach  is  to  modify  a  surface  near  the  source  of  the  blast 
waves  to  reduce  the  amount  of  blast  wave  energy  propagating  intothe  far  field.  By 
using  a  source  above  an  absorbing  surface,  the  energy  reflected  from  the  ground  is 
reduced.  Because  a  rigid  ground  surface  enhances  the  sound  level  by  6  dB,  reduc¬ 
ing  or  eliminating  this  energy  can  potentially  have  a  significant  impact  on  sound 
levels.  Also,  by  shielding  a  source  with  an  absorbing  surface,  blast  wave  energy  is 
absorbed. 

A  previous  ERDC  research  project  showed  that  a  gravel  pit  beneath  the  source 
location  could  reduce  the  sound  levels  produced  by  small  explosions  (Attenbor¬ 
ough  et  al.  2003).  However,  there  were  some  inconsistencies  in  these  tests,  and 
the  tests  did  not  include  measurements  at  longer  distances  that  would  be  important 
in  practice.  In  addition,  a  better  design  of  the  blast  absorber  might  provide  addi¬ 
tional  attenuation.  This  field  test  was  designed  to  answer  these  remaining  ques¬ 
tions. 

The  previous  project  was  concerned  primarily  with  developing  theoretical  and 
computational  methods  to  understand  and  predict  blast  absorption  from  granular 
materials  at  high  acoustic  pressures.  The  absorption  pit  sizes  were  limited  by  cost, 
and  were  not  optimized  to  reduce  low  frequency  blast  noise.  The  gravel  pit  used 
in  the  earlier  measurements  was  3.66  x  3.66  m,  1.5  m  deep  (12  x  12  x  5  ft  deep) 
and  filled  with  either  three-layered  gravel  or  a  uniform  gravel  layer.  The  pit  with  a 
uniform  0.9-cm  gravel  size  seemed  to  produce  a  greater  reduction  in  blast  noise 
that  the  pit  with  three  layers. 

To  improve  the  absorber  performance  at  low  frequencies,  the  size  of  the  gravel  pit 
used  for  blast  absorption  must  be  increased  to  roughly  a  wavelength  in  size.  Be¬ 
cause  a  single  brick  charge  of  C-4  typically  has  a  peak  frequency  of  30-50  Hz,  the 
wavelength  is  roughly  8-12  m.  Thus,  this  study  chose  a  50-  x  50-ft  (approxi¬ 
mately  15-  x  15-m)  size  for  the  gravel  pit.  Because  of  the  presence  of  unexploded 
ordinance  on  the  test  range,  it  was  not  possible  to  dig  a  pit  and  fill  it  with  gravel. 
Instead  a  pile  of  gravel,  1.5-m  high  with  steep  edges  of  approximately  45°,  was 
constructed  instead  of  creating  a  pit  in  the  ground.  As  the  previous  measurements 
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showed  a  large  reduction  (40  dB)  in  the  blast  wave  after  penetration  to  a  depth  of 
0.5  m,  we  also  specified  a  depth  of  5  ft  (1.5  m)  for  the  pad  as  it  will  be  enough  to 
significantly  reduce  any  reflections  from  the  bottom  of  the  pad.  Coarse  gravel  (3 
cm)  was  used  to  create  the  pad  and  a  variety  of  source  locations  were  used  for 
comparison  (Figures  1  and  2). 


Figure  1.  The  1.5-m-height  gravel  pad. 
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3  Instrumentation  Sensor  Array 

Two  different  sensor  arrays  were  used.  To  determine  the  effect  of  the  blast  ab¬ 
sorber  at  relatively  short  propagation  distances,  ERDC-CRREL  installed  an  array 
of  36  sensors  located  from  12  to  400  m  away  from  the  source  locations  (Figure  3). 
The  origin  (reference  point)  is  even  with  the  back  side  of  the  gravel  pad,  in  line 
with  the  gravel  pad  array.  Because  of  the  test  geometry,  two  different  pressure 
sensors  were  installed  to  record  the  noise  levels  at  the  closest  distances,  one  for 
each  source  location.  For  distances  greater  than  100  m,  the  same  sensors  were 
used  at  each  distance.  The  closest  sensor  locations  at  12,  22,  and  32  m  have  pres¬ 
sure  sensors  only,  as  the  ground  motion  sensors  would  be  overdriven  and  would 
not  produce  valid  data.  Beyond  32  m,  each  sensor  location  had  two  pressure  sen¬ 
sors  (microphones),  one  at  the  ground  surface  and  one  a  1.5-m  height,  and  two 
geophones  to  measure  the  vertical  and  horizontal  ground  motions.  The  linear  sen¬ 
sor  array  provided  measurements  at  distances  of  32,  62,  92,  120,  180,  240,  340, 
370  and  400  m  from  each  source  location. 


The  second  sensor  array  (ERDC-CERL)  was  used  to  determine  the  blast  noise  at 
longer,  more  realistic  distances  of  1.8,  3.2  and  4.5  km.  These  measurements  will 
be  discussed  in  a  separate  report. 
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Figure  3.  ERDC-CRREL  instrumentation  array  and  C-4  detonation  locations. 


All  the  shot  points  and  sensor  stations  locations  are  presented  in  Table  1 .  All 
measurements  are  in  meters.  Table  2  gives  information  about  the  ERDC-CRREL 
sensor  array,  while  Figure  4  shows  a  typical  instrumentation  station. 
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Table  1.  Survey  coordinates  of  shot  points  and  sensor  stations. 


TW  #1  Blast  Absorber  Test 
Coordinate  System:  UTM 


Datum:  NAD83 
Units:  Meters 

North 

East 

Description 

Grid  N 

Grid  E 

Distance 

2  mL 

2  mR 

-5  mL 

-5  mR 

4368677.645 

402583.301 

"0"  REF 

0.0 

0.0 

0.0 

15.2 

15.1 

15.8 

15.8 

4368673.932 

402598.677 

-5  m  bare  gnd. 

-3.7 

15.4 

15.8 

7.0 

30.8 

0.0 

30.0 

4368668.001 

402594.991 

2  m  bare  gnd. 

-9.6 

11.7 

15.2 

0.0 

30.0 

7.0 

30.8 

4368663.758 

402592.347 

7  m  bare  gnd 

-13.9 

9.0 

16.6 

5.0 

30.4 

12.0 

32.3 

4368659.525 

402589.700 

12  m  bare  gnd. 

-18.1 

6.4 

19.2 

10.0 

31.6 

17.0 

34.5 

4368651.036 

402584.395 

22  m  bare  gnd. 

-26.6 

1.1 

26.6 

20.0 

36.1 

27.0 

40.4 

4368642.556 

402579.094 

32  m  bare  gnd. 

-35.1 

-4.2 

35.3 

30.0 

42.4 

37.0 

47.6 

4368617.078 

402563.229 

62  m  bare  gnd. 

-60.6 

-20.1 

63.8 

60.0 

67.1 

67.0 

73.5 

4368591.630 

402547.348 

92  m  bare  gnd. 

-86.0 

-36.0 

93.2 

90.0 

94.9 

97.0 

101.6 

4368607.501 

402521.883 

92  m  gravel  pad 

-70.1 

-61.4 

93.2 

94.9 

90.0 

101.5 

97.0 

4368632.952 

402537.755 

62  m  gravel  pad 

-44.7 

-45.5 

63.8 

67.1 

60.0 

73.4 

67.1 

4368658.424 

402553.640 

32  m  gravel  pad 

-19.2 

-29.7 

35.3 

42.4 

30.0 

47.6 

37.0 

4368666.898 

402558.948 

22  m  gravel  pad 

-10.7 

-24.4 

26.6 

36.1 

20.0 

40.3 

27.0 

4368675.408 

402564.232 

12  m  gravel  pad 

-2.2 

-19.1 

19.2 

31.6 

10.0 

34.5 

17.0 

4368679.641 

402566.890 

7  m  gravel  pad 

2.0 

-16.4 

16.5 

30.4 

5.0 

32.3 

12.0 

4368683.893 

402569.523 

2  m  gravel  pad 

6.2 

-13.8 

15.1 

30.0 

0.0 

30.8 

7.0 

4368689.843 

402573.246 

-5  m  gravel  pad 

12.2 

-10.1 

15.8 

30.8 

7.0 

30.0 

0.0 

4368569.529 

402515.876 

120  m 

-108.1 

-67.4 

127.4 

126.3 

126.3 

133.3 

133.3 

4368524.931 

402488.029 

180  m 

-152.7 

-95.3 

180.0 

178.6 

178.6 

185.6 

185.6 

4368467.713 

402452.364 

240  m 

-209.9 

-130.9 

247.4 

245.9 

245.9 

252.9 

252.9 

4368382.869 

402399.438 

340  m 

-294.8 

-183.9 

347.4 

345.7 

345.8 

352.7 

352.8 

4368363.386 

402387.299 

370  m 

-314.3 

-196.0 

370.4 

368.7 

368.7 

375.7 

375.7 

4368337.850 

402371.359 

400  m 

-339.8 

-211.9 

400.5 

398.8 

398.8 

405.7 

405.8 

4368252.771 

402318.232 

500  m 

-424.9 

-265.1 

500.8 

499.0 

499.0 

506.0 

506.0 

4367826.416 

402052.341 

1000  m 

-851.2 

-531.0 

1003.2 

1001.4 

1001.4 

1008.3 

1008.4 

4367152.301 

401536.859 

1845  m 

-1525.3 

-1046.4 

1849.8 

1848.5 

1847.2 

1855.5 

1854.2 

4365918.666 

400925.926 

3218.5  m 

-2759.0 

-1657.4 

3218.5 

3216.3 

3216.8 

3223.3 

3223.8 

4364792.940 

400394.874 

4458.7  m 

-3884.7 

-2188.4 

4458.7 

4456.1 

4457.4 

4463.1 

4464.4 

Table  2.  ERDC-CRREL  sensor  array. 


a.  Cable  #1,  30-m  array  cable  (four  takeouts  every  30  m),  from  gravel  pad  source  location  (at 
+2  m). 


Station 

Channel 

Sensor 

S/N 

Location  (m) 

Height  (cm) 

1 

12 

PCB  102A06 

15962 

12 

10.2 

1 

11 

PCB  102A06 

15964 

22 

152.4 

1 

10 

PCB  102A09 

15969 

32 

152.4 

1 

9 

PCB  102A09 

15968 

32 

12.7 

2 

8 

PCB  102A07 

15971 

62 

152.4 

2 

7 

PCB  102A07 

15973 

62 

16.5 

2 

6 

V  geophone 

62 

0 

2 

5 

H  geophone 

62 

0 

3 

3 

PCB  106B50 

6522 

92 

152.4 

3 

4 

PCB  106B50 

6693 

92 

10.8 

3 

2 

V  geophone 

92 

0 

3 

1 

H  geophone 

92 

0 
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Table  2  (cont’d).  ERDC-CRREL  sensor  array. 

b.  Cable  #2,  30-m  array  cable  (four  takeouts  every  30  m),  from  bare  ground  source  location. 


Station 

Channel 

Sensor 

S/N 

Location  (m) 

Height  (cm) 

4 

24 

PCB  102A06 

15965 

12 

8.9 

4 

23 

PCB  102A06 

15966 

22 

142.2 

4 

22 

PCB  102A09 

15970 

32 

148.6 

4 

21 

PCB  102A09 

13180 

32 

15.2 

5 

20 

PCB  102A07 

15972 

62 

151.1 

5 

19 

PCB  102A07 

15974 

62 

9.5 

5 

18 

V  geophone 

62 

0 

5 

17 

H  geophone 

62 

0 

6 

16 

PCB  106B50 

6695 

92 

152.4 

6 

15 

PCB  106B50 

3259 

92 

16.5 

6 

14 

V  geophone 

92 

0 

6 

13 

H  geophone 

92 

0 

c.  Cable  #3,  60-m  interval  between  every  four  takeouts.  Set  attenuation  factor  to  1/10  (-20 
dB)  for  all  B&K  microphones  on  this  cable. 


Station 

Channel 

Sensor 

S/N 

Location  (m) 

Height  (cm) 

7 

36 

B&K  4938-1 

2239253 

120 

152.4 

7 

35 

B&K  4938-2 

2239254 

120 

25.4 

7 

34 

V  geophone 

120 

0 

7 

33 

H  geophone 

120 

0 

8 

32 

B&K  4938-3 

2239255 

180 

152.4 

8 

31 

B&K  4938-4 

2239256 

180 

25.4 

8 

30 

V  geophone 

180 

0 

8 

29 

H  geophone 

180 

0 

9 

28 

B&K  4938-5 

2239257 

240 

150 

9 

27 

B&K  4938-6 

2239258 

240 

25.4 

9 

26 

V  geophone 

240 

0 

9 

25 

H  geophone 

240 

0 
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Table  2  (cont’d).  ERDC-CRREL  sensor  array. 


d.  Cable  #4,  200-m  interval  between  every  four  takeouts.  Set  attenuation  factor  to  1/100  (- 
40  dB)  for  all  B&K  microphones  on  this  cable. 


Station 

Channel 

Sensor 

S/N 

Location  (m) 

Height  (cm) 

10 

37 

B&K  4165-7 

1881046 

340 

150 

10 

47 

B&K  4165-8 

1881049 

340 

0 

10 

46 

V  geophone 

340 

0 

10 

45 

H  geophone 

340 

0 

11 

44 

B&K  4165-9 

1881050 

370 

150 

11 

43 

B&K  4165-10 

2068999 

370 

0 

11 

42 

V  geophone 

370 

0 

11 

41 

H  geophone 

370 

0 

12 

40 

B&K  4165-5 

1857589 

400 

150 

12 

39 

B&K  4165-6 

1881043 

400 

0 

12 

38 

V  geophone 

400 

0 

12 

37 

H  geophone 

400 

0 

Microphone  at  1.5m  lieieht 


Microphone  at  ground  surface 


Vertical  and  radial  gecphone* 


Figure  4.  Instrumentation  Station  6  at  90  m. 
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4  Recorded  Events 

Data  was  recorded  on  a  digital  Strata  Visor  (Figure  5  and  Table  3)  NZ  Exploration 
Seismic  Recorder  (Appendix  B),  with  48  channels  and  0.03125  ms  sampling  in¬ 
terval. 


Figure  5.  StrataVisor  NZ  seismograph  recorder. 


Table  3.  StrataVisor  Geometric  Model  NZ  seismograph  specifications. 


Number  of  channels 

48 

Digitization  (bits) 

24 

Sampling  rate 

0.03125  ms 

Frequency  band 

1.75  Hz  -  20  kHz  (-3dB  at  83%  Nyquist,  Down  90  dB 
at  Nyquist) 

Maximum  signal  voltage 

±1.4  V  at  0  dB 

Maximum  record  length 

65k 

Noise  floor,  jliV  RMS 

0.2 

Pretrigger 

4096  samples 

Operating  temperature 

+5°C  start-up,  -5°C  operation 

Records  were  obtained  for  thirty- three  C-4  detonations,  twenty  .45 -caliber  blank 
pistol  shots  and  six  detonated  blasting  caps.  Numerous  noise  records  and  micro¬ 
phone  calibration  records  were  also  taken.  Only  the  C-4,  pistol,  and  blasting  cap 
records  are  presented  in  this  report. 
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C-4  Recorded  Events. 

Thirty-three  single  brick  Ml  12  C-4  charges  were  detonated  at  six  different  shot 
points.  The  bricks  weigh  1.25  lb  (0.57  kg)  and  were  detonated  above  the  ground 
by  suspending  them  from  a  cross  line  between  two  poles  (Figure  6). 


wHWsi.V  * 

Figure  6.  C-4  brick  being  suspended  at  1.5  m  above  the  gravel  pad. 


The  different  shot  point  locations  were  identified  by  a  code  of  the  form  Type  + 
Distance  +  Height.  Here,  Type  is  a  letter  code  denoting  whether  the  shot  was  near 
the  absorber  or  not.  P  is  the  letter  used  to  specify  the  gravel  pad,  while  B  indicates 
that  the  shot  was  over  bare  ground.  The  Distance  is  given  by  a  number,  either  +2 
or  -5.  This  is  the  x  grid  coordinate  in  meters,  with  0  being  the  back  edge  of  the 
gravel  pad.  Finally,  H  is  the  height  of  the  source  above  the  ground  or  gravel  pad  in 
meters  (Figure  7).  Table  4  lists  the  events. 


Figure  7.  Locations  of  the  six  C-4  shot  points. 
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Table  4.  C-4  recorded  events  (shot  records). 


CRREL 
Binary 
File  No. 

CERL 

Shot 

No. 

Date 

Time 

(EST) 

Source 

Source 

Location 

Source  Remarks 

Height 

(m) 

37 

1 

15-Jun 

13:46 

C4 

B-5+3 

Ch  21  &  15  bad,  18,  11  &  12  have  crosstalk 

38 

2 

15-Jun 

14:42 

C4 

B+2+1.5 

39 

3 

15-Jun 

14:34 

C4 

P+2+1.5 

(FP2  back)  Ch  1 1  &  12  timebreak  crosstalk,  Ch  21  &  15  bad 

40 

4 

15-Jun 

15:55 

C4 

B+2+1.5 

(FP1  back)  Ch  21  bad 

41 

15-Jun 

Noise  record 

42 

5 

15-Jun 

15:58 

C4 

P+2+1.5 

(FP2  back)  Ch  21  bad.  Swapped  out  Ch21 

W/102A06,  s/n  16784,  30m 

43 

6 

15-Jun 

16:31 

C4 

B-5+3 

Good  record,  all  channels 

44 

7 

15-Jun 

16:33 

C4 

P-5+3 

Good  record,  all  channels  except  11  &  12 

45 

8 

15-Jun 

16:55 

C4 

B-5+3 

Cross  talk  on  16 

46 

9 

15-Jun 

16:57 

C4 

P-5+3 

Cross  talk  on  18,  12  looks  noisy 

68 

16-Jun 

9:39 

Noise  record 

69 

10 

16-Jun 

9:50 

C4 

B-5+1 

Ch  18  &  15  early  crosstalk 

70 

11 

16-Jun 

9:52 

C4 

P-5+1 

Ch  18,  15  &  23  crosstalk 

71 

12 

16-Jun 

9:45 

C4 

B-5+1 

Ch  18,  15  &  23  crosstalk 

72 

13 

16-Jun 

10:11 

C4 

P-5+1 

Ch  18  &  15  crosstalk 

16-Jun 

10:15 

Moved  cable  near  Ch  18 

73 

14 

16-Jun 

10:30 

C4 

B+5+1.5 

Ch  18  &  15  crosstalk 

74 

15 

16-Jun 

10:31 

C4 

P-5+1.5 

Ch  18  &  15  crosstalk 

75 

16 

16-Jun 

10:46 

C4 

B+2+1.5 

Ch  18  crosstalk 

76 

17 

16-Jun 

10:48 

C4 

P+2+1.5 

Ch  18  &  15  crosstalk 

77 

18 

16-Jun 

11:01 

C4 

B+2+1.5 

Ch  18  crosstalk,  a  little  on  15 

78 

19 

16-Jun 

11:03 

C4 

P+2+1.5 

Ch  18  &  15  crosstalk 

79 

20 

16-Jun 

11:18 

C4 

B+2+1.5 

Ch  18  &  15  crosstalk 

80 

21 

16-Jun 

11:19 

C4 

P+2+1.5 

Ch  18,  15  &  23  crosstalk 

81 

22 

16-Jun 

11:43 

C4 

B-5+3 

Ch  18  crosstalk 

82 

23 

16-Jun 

11:45 

C4 

P-5+3 

AQ  little  crosstalk  on  Ch  18 

83 

24 

16-Jun 

13:53 

C4 

B-5+3 

Good  record 

84 

25 

16-Jun 

13:54 

C4 

P-5+3 

A  little  noise  on  Ch  23 

85 

26 

16-Jun 

14:08 

C4 

B-5+3 

All  channels  good 

86 

27 

16-Jun 

14:09 

C4 

P-5+3 

All  channels  good 

87 

28 

16-Jun 

14:25 

C4 

B-5+1 

TD  crosstalk  on  many  channels 

88 

29 

16-Jun 

14:26 

C4 

P-5+1 

TD  crosstalk  on  many  channels 

89 

30 

16-Jun 

14:40 

C4 

B-5+1 

Some  crosstalk  on  CH  11,  12.  15,  21  &  23 

90 

31 

16-Jun 

14:42 

C4 

P-5+1 

Some  crosstalk  on  CH  15  &  23 

91 

32 

16-Jun 

14:55 

C4 

B-5+1 

TD  crosstalk  on  many  channels 

92 

33 

16-Jun 

14:56 

C4 

P-5+1 

Some  crosstalk  on  multiple  channels 

Blank  Pistol  Recorded  Events 

A  .45-caliber  blank  pistol  was  held  at  a  height  of  1  m  and  fired  along  the  array  to 
characterize  the  ground  impedance  (Figure  8).  The  blank  cartridges  used  are  cus¬ 
tom  made  with  an  overall  length  of  22.6  mm,  7.6-g  cartridge  weight,  5.8-g  case 
weight,  with  31  g  of  FFFFG,  CCI  black  powder.  Figure  9  shows  location  and  di¬ 
rection  of  the  blank  pistol  shots,  while  Table  5  lists  the  events. 
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Figure  8.  Blank  .45-caliber  pistol  fired  along  the  sensor  array  at  a  1-m  height. 


Figure  9.  Location  and  direction  of  the  blank  pistol  shots. 
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Table  5.  Blank  pistol  recorded  events  (shot  record). 


CRREL 
Binary 
File  No. 

CERL 

Shot 

No. 

Date 

Time 

(EST) 

Source 

Source 

Location 

Source 

Height 

(m) 

Remarks 

47 

16-Jun 

8:19 

Noise  record 

NZ  @.  125ms,  pretrigger -.5sec 

48 

16-Jun 

8:23 

Pistol  1 

280m,  up  line 

1 

Ch  47  &  48  were  off 

49 

16-Jun 

Pistol  1 

280m,  up  line 

1 

50 

16-Jun 

Pistol  1 

280m,  up  line 

1 

51 

16-Jun 

Pistol  1 

280m,  up  line 

1 

52 

16-Jun 

8:32 

Pistol  1 

280m,  up  line 

1 

53 

16-Jun 

Pistol  2 

310m,  up  line 

1 

54 

16-Jun 

Pistol  2 

310m,  up  line 

1 

55 

16-Jun 

Pistol  2 

310m,  up  line 

1 

56 

16-Jun 

Pistol  2 

310m,  up  line 

1 

57 

16-Jun 

Pistol  2 

310m,  up  line 

1 

58 

16-Jun 

Pistol  3 

430m,  down  line 

1 

59 

16-Jun 

8:45 

Pistol  3 

430m,  down  line 

1 

60 

16-Jun 

Pistol  3 

430m,  down  line 

1 

61 

16-Jun 

Pistol  3 

430m,  down  line 

1 

62 

16-Jun 

8:48 

Pistol  3 

430m,  down  line 

1 

63 

16-Jun 

8:53 

Pistol  4 

280m,  down  line 

4 

Personnel  @  90m 

64 

16-Jun 

Pistol  4 

280m,  down  line 

4 

65 

16-Jun 

Pistol  4 

280m,  down  line 

4 

66 

16-Jun 

Pistol  4 

280m,  down  line 

4 

67 

16-Jun 

8:57 

Pistol  4 

280m,  down  line 

4 

Blasting  Cap  Recorded  Events 

Six  blasting  caps  were  detonated  from  two  different  shot  points  to  assist  in  char¬ 
acterizing  the  ground  impedance  and  instrumentation  response.  The  caps  were 
suspended  above  the  ground  in  a  manner  similar  to  the  C-4  charges.  Figure  10 
gives  the  locations  of  the  blasting  cap  detonations,  while  Table  6  lists  the  events. 


Table  6.  Blasting  cap  recorded  events  (shot  record). 


CRREL 
Binary 
File  No. 

CERL 

Shot 

No. 

Date 

Time 

(EST) 

Source 

Source 

Location 

Source  Remarks 

Height 

(m) 

105 

16-Jun 

15:38 

Noise  record 

106 

Cap  1 

16-Jun 

15:44 

Blasting  cap 

P+2+2 

Switched  to  36  channels  on  NZ 

107 

Cap  2 

16-Jun 

15:49 

Blasting  cap 

P+2+2 

108 

Cap  3 

16-Jun 

16:01 

Blasting  cap 

B+2+2 

Crosstalk 

109 

Cap  4 

16-Jun 

16:25 

Blasting  cap 

B+2+2 

Crosstalk 

110 

Cap  5 

16-Jun 

16:32 

Blasting  cap 

B+2+2 

111 

Cap  6 

16-Jun 

16:37 

Blasting  cap 

B+2+2 

Lots  of  crosstalk 
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5  Meteorological  Measurements  During  Testing. 

The  ERDC-CRREL  met  station  (Figure  11),  a  Hobo  model  H-08-04-002,  was  lo¬ 
cated  approximately  10  m  off  set  from  Station  7  at  120  m.  Measurements  of  tem¬ 
perature,  wind  speed,  wind  gusts,  wind  direction,  atmospheric  pressure,  and  rela¬ 
tive  humidity  were  recorded  during  the  three  days  of  testing.  Figures  12  and  13 
show  the  meteorological  data  during  the  shot  testing  times  on  16  and  17  June. 
Figure  14  shows  the  complete  3 -day  meteorological  data  over  the  entire  time  dur¬ 
ing  shots  and  calibrations. 


Figure  11.  ERDC-CRREL  Met  Station  at  120  m. 


Pressure  (mbar) 
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Pressure  (mbar) 
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Figure  14.  Met  data  for  the  entire  3-day  testing  period. 
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6  Data  Quality,  Signature  Records,  and  Data 
Plots. 

Each  channel  from  all  the  records  is  evaluated  for  response,  noise,  ringing,  and 
crosstalk  to  determine  if  that  data  is  usable.  The  following  three  Tables  7  through 
9  rate  each  data  channel. 


Table  7.  Pistol  sensor  data  quality. 


MD  2005  Blast  Absorber  Test 

Blank  Pistol  Shots 

Sensors  -  Channel  numbers 

Note:  Channels  1-24  are  bad  for  all  pistol  shots 

CRREL  Rec  Num-> 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67  68  69 

36 

P-122-1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

VN 

X 

N 

N 

O 

O 

O 

35 

P-122-0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

VN 

X 

N 

N 

N 

VN 

N 

34 

V-122 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

VN 

X 

N 

33 

R-122 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

VN 

VN 

N 

X 

N 

32 

P-182-1 

VN 

VN 

VN 

VN 

VN 

X 

X 

X 

X 

X 

X 

VN 

VN 

VN 

VN 

O 

O 

O 

o 

O 

31 

P-182-0 

VN 

VN 

VN 

VN 

VN 

X 

X 

X 

X 

X 

X 

VN 

VN 

VN 

VN 

o 

o 

O 

o 

O 

30 

V-182 

VN 

VN 

VN 

VN 

VN 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

N 

N 

o 

N 

o 

29 

R-182 

VN 

VN 

VN 

VN 

VN 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

N 

N 

o 

N 

o 

28 

P-242-1 

O 

O 

o 

o 

o 

VN 

VN 

VN 

VN 

VN 

VN 

VN 

VN 

VN 

VN 

O 

O 

o 

o 

o 

27 

P-242-0 

o 

o 

o 

o 

o 

VN 

VN 

VN 

VN 

VN 

VN 

VN 

VN 

VN 

VN 

o 

o 

o 

o 

o 

26 

V-242 

VN 

VN 

VN 

VN 

VN 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

o 

o 

o 

o 

o 

25 

R-242 

N 

N 

N 

N 

N 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

o 

o 

o 

o 

o 

48 

P-342-1 

X 

O 

O 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

VN 

47 

P-342-0 

X 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

N 

46 

V-342 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

N 

o 

o 

N 

N 

N 

N 

N 

45 

R-342 

o 

N 

N 

N 

N 

o 

o 

o 

o 

o 

N 

o 

N 

N 

N 

VN 

VN 

VN 

VN 

VN 

44 

P-372-1 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

N 

N 

N 

N 

N 

43 

P-372-0 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

N 

N 

N 

N 

N 

42 

V-372 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

N 

o 

o 

N 

N 

N 

VN 

VN 

41 

R-372 

o 

o 

o 

o 

o 

o 

o 

o 

N 

o 

N 

o 

N 

N 

N 

N 

X 

VN 

VN 

VN 

40 

P-402-1 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

VN 

VN 

VN 

VN 

N 

39 

P-402-0 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

VN 

VN 

VN 

VN 

N 

38 

V-402 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

VN 

VN 

VN 

VN 

VN 

37 

R-402 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

VN 

X 

VN 

VN 

VN 

1-12 

ALL 

CHANNELS  BAD  ON  ALL  PISTOL  SHOTS 

13-24 

ALL  CHANNELS  BAD  ON  ALL  PISTOL  SHOTS 

O  -  Good 

X  =  Bad 

N  =  Noisy  data  -  but  usable 

R  =  Ringing  sensor 

C  =  Crosstalk 

VN  =  Very  Noisy  data  - 

arrival  visible  but  not  usable  without  add'l  processing 
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Table  8.  C-4  sensor  data  quality. 


MD  2005  Blast  Absorber  Test 

C-4  Shots 

Sensors  -  Channel  numbers 

CERL  Shot  Num  ■> 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

CRREL  Rec  Num-> 

37 

38 

39 

40 

42 

43 

44 

45 

46 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

Source  Location  ■> 

B 

B 

p 

B 

P 

B 

P 

B 

P 

B 

p 

B 

p 

B 

P 

B 

P 

B 

P 

B 

P 

B 

P 

B 

P 

B 

P 

B 

P 

B 

P 

B 

P 

12 

PP-12-0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

O 

0 

0 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

11 

PP-22-1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

O 

0 

0 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

O 

O 

O 

O 

O 

O 

0 

O 

0 

O 

0 

10 

PP-32-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

9 

PP-32-0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

8 

PP-62-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

7 

PP-62-0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

6 

VP-62 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

5 

RP-62 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

4 

PP-92-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

O 

0 

O 

0 

N 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

3 

PP-92-0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

O 

0 

2 

VP-92 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

0 

0 

0 

0 

0 

0 

O 

0 

0 

0 

O 

0 

1 

RP-92 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

0 

0 

O 

0 

24 

PB-12-0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

23 

PB-22-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

22 

PB-32-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

21 

PB-32-0 

X 

X 

X 

X 

X 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

20 

PB-62-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

19 

PB-62-0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

N 

0 

0 

0 

0 

0 

0 

18 

VB-62 

X 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

17 

RB-62 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

16 

PB-92-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

15 

PB-92-0 

X 

X 

X 

0 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c 

c 

c 

c 

c 

0 

N 

0 

0 

N 

N 

N 

0 

14 

VB-92 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

13 

RB-92 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

36 

P-122-1 

0 

35 

P-122-0 

0 

34 

V-122 

0 

33 

R-122 

0 

32 

P-182-1 

0 

ALL  CHANNELS  GOOD  ON  ALL  SHOTS 

31 

P-182-0 

0 

30 

V-182 

0 

29 

R-182 

0 

28 

P-242-1 

0 

27 

P-242-0 

0 

26 

V-242 

0 

25 

R-242 

0 

48 

P-342-1 

0 

47 

P-342-0 

0 

46 

V-342 

0 

45 

R-342 

0 

44 

P-372-1 

0 

ALL  CHANNELS  GOOD  ON  ALL  SHOTS 

43 

P-372-0 

0 

42 

V-372 

0 

41 

R-372 

0 

40 

P-402-1 

0 

39 

P-402-0 

0 

38 

V-402 

0 

37 

R-402 

0 

O  =  Good 

X  =  Bad 

N  =  Noisy  data  -  but  usable 

R  =  Ringing  sensor 

C  =  Crosstalk 
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Table  9.  Blasting  cap  sensor  data  quality. 

MD  2005  Blast  Absorber  Test  Blasting  caps 

Sensors  -  Channel  numbers 


CRREL  Rec  Num-> 

106 

107 

108 

109 

110 

111 

12 

PP-12-0 

O 

O 

O 

X 

X 

X 

11 

PP-22-1 

VN 

N 

N 

X 

X 

X 

10 

PP-32-1 

O 

O 

O 

VN 

VN 

o 

9 

PP-32-0 

VN 

O 

O 

VN 

VN 

N 

8 

PP-62-1 

VN 

o 

o 

X 

X 

X 

7 

PP-62-0 

O 

o 

o 

VN 

VN 

VN 

6 

VP-62 

o 

o 

o 

O 

O 

N 

5 

RP-62 

o 

o 

o 

o 

o 

N 

4 

PP-92-0 

VN 

N 

o 

N 

VN 

X 

3 

PP-92-1 

N 

o 

o 

X 

X 

X 

2 

VP-92 

o 

o 

o 

o 

o 

N 

1 

RP-92 

o 

o 

o 

o 

o 

N 

24 

PB-12-0 

X 

X 

VN 

o 

N 

N 

23 

PB-22-1 

X 

X 

N 

N 

VN 

N 

22 

PB-32-1 

X 

X 

VN 

N 

N 

N 

21 

PB-32-0 

X 

X 

X 

X 

X 

VN 

20 

PB-62-1 

X 

X 

VN 

VN 

X 

X 

19 

PB-62-0 

X 

X 

VN 

X 

VN 

X 

18 

VB-62 

o 

o 

O 

o 

o 

N 

17 

RB-62 

o 

o 

o 

o 

o 

O 

16 

PB-92-1 

o 

N 

o 

VN 

VN 

VN 

15 

PB-92-0 

VN 

N 

N 

VN 

VN 

X 

14 

VB-92 

o 

N 

o 

o 

o 

N 

13 

RB-92 

o 

o 

o 

o 

o 

N 

36 

P-122-1 

N 

N 

N 

N 

N 

VN 

35 

P-122-0 

N 

N 

N 

N 

N 

VN 

34 

V-122 

N 

O 

O 

O 

N 

N 

33 

R-122 

N 

N 

o 

o 

N 

X 

32 

P-182-1 

VN 

O 

o 

o 

X 

N 

31 

P-182-0 

VN 

O 

o 

o 

X 

N 

30 

V-182 

O 

o 

o 

o 

X 

O 

29 

R-182 

o 

o 

o 

o 

X 

O 

28 

P-242-1 

X 

VN 

N 

VN 

VN 

X 

27 

P-242-0 

X 

VN 

N 

N 

VN 

X 

26 

V-242 

X 

VN 

N 

N 

VN 

X 

25 

R-242 

X 

VN 

N 

N 

VN 

X 

O  =  Good 
X  =  Bad 

N  =  Noisy  data  -  but  usable 
R  =  Ringing  sensor 
C  =  Crosstalk 

VN  =  Very  Noisy  data  -  arrival  visible  but  not  usable  without  add’l  processing 


Note:  Channels  36-48  were  off  during  these  recordings  (except  record  106,  bad) 

108:  Cross-talk  interference  on  all  channels  -  muting  needed 
112:  Strong  time  break  cross-talk  on  all  channels 

The  following  signal  data  plots  show  the  response  of  the  sensors  recorded  (48  for 
C-4  and  the  blasting  caps  and  24  for  the  pistol)  during  each  test.  The  plots  are 
normalized  so  that  each  channel  is  the  same  size.  The  maximum  amplitude  (in  Pa 
for  pressure  sensors,  m/s  for  geophones)  is  shown  in  the  label  for  each  channel. 
Figures  15  through  17  describe  the  plot  layout,  labels,  and  descriptions  for  the  C- 
4,  pistol,  and  blasting  cap  plots. 
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Figure  15.  Description  of  C-4  data  plots. 
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!  Project  Name 


Record  Number 


Source 


Source  Location 


'4 


MD  2005  —  Record  48  Pistol  280  m 


160-0.43 
160-0.36 
160-1.2E-6 
160-  3.  IE-6 
100-0.88 
100-0.63 
100-  6.2E-6 
100-  5.7E-6 
40-4.39 
40-1.87 
40-  3.3E-6 
40-  8.6E-6 
60-0.45 
60-97.7E-3 
60-18.2E-6 
60-13.0E-6 
90-5.92 
90-1.86 
90-12.7E-6 
90-15.6E-6 
120-1.59 
120-1.21 
120-  9.5E-6 
120-  7.9E-6 


P-120-1 


P-120-0 


R-120 


P-180-0 
V-180 
R-180 


4- 


Ii^w l 


jut 


***#**0#*S  V-240 

"**"*~~~~»*  R-240 


I** — . - . — — — — — 


-p 


~||jUw 


LL— .. 


P-180-1 


P-240-1 


P-240-0 


P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 

t 


*  -0.5 


Location  &  Peak  amplitude 

(Pa  for  pressure) 

(mm/s  for  geophones) 


0.5 


1.5 


Time,  sec 


2f 


Sensor  Name  R  -  400 ' 


Sensor  type 

P  =  pressure 
V  =  vertical  geophone 
R  =  radial  geophone 


Sensor  location 

Other  examples: 

10,  30,  60,  90  up  to 
400m  (from  source) 


Figure  16.  Description  of  pistol  data  plots. 
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MD  2005  —  Record  106  Cap  P+2+2 


10-963.0 
20-304.0 
30-248.0 
30-72.4 
60-105.8 
60-86.3 
60-80.5E-6 
60-96.3E-6 
90-10.5 
90-26.0 
90-27. OE-6 
90-60. OE-6 
32-358.8 
36-806.7 
42-119.5 
42-157.7 
67-73.5 
67-92.6 
67-59.7E-6 
67-89.2E-6 
95-74.5 
95-27.8 
95-1 8. OE-6 
95-27.7E-6 
119-8.60 
119-5.73 
119-21.6E-6 
1 19-32.5E-6 
179-3.54 
179-3.09 
179-30.  IE-6 
1 79-56. 2E-6 
238-2.46 
238-2.10 
238-  6. IE-6 
238-21. 4E-6 
338-3.01 
338-0.71 
338-  6.7E-6 
338-1 1.6E-6 
368-3.05 
368-1 .03 
368-12.  IE-6 
368- 13.  IE-6 
398-6.41 
398-2.64 
398-  9.2E-6 
398- 13.  IE-6 


H— h 


PP-10-0 

PP-22-1 

PP-32-1 

PP-32-0 

PP-62-1 

PP-62-0 

VP-62 

RP-62 

PP-92-0 

PP-92-1 

VP-92 

RP-92 

PB-10-0 

PB-22-1 

PB-32-1 

PB-32-0 

PB—62-1 

PB-62-0 

VB-62 

RB-62 

PB-92-1 

PB-92-0 

VB-92 

RB-92 

P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 


H  0  02  0.4  0.6 

* 

% 

% 

0.8  1  1.2  1.4 

Time,  sec 

1.6  1.8  f 

t 

1 

1 

J"  i 

,  Location  &  Peak  amplitude  | 

Sensor  Name  R 

-  400 

!  (Pa  for  pressure)  J 

!  (mm/s  for  geophones)  [ 

!  1 

Sensor  type 

Sensor  location 

P  =  pressure 

Other  examples: 

V  =  vertical  geophone 

10 ,  30?  60, 90  up  to 

R  =  radial  geophone 

400m  (from  source) 

Figure  17.  Description  of  blasting  cap  data  plots. 
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Pistol  Data  Plots 

MD  2005  — Record  48  Pistol  280  m 

P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 
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V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 

-0.5  0  0.5  1  1.5  2 

Time,  sec 


160-0.43 

160-0.36 

160-  1.2E-6 

160- 3. IE-6 

100-0.88 

100-0.63 
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40-4.39 

40-1.87 

40-  3.3E-6 

40-  8.6E-6 

60-0.45 
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60-18.2E-6 

60-13. OE-6 

90-5.92 

90-1.86 

90-12. 7E-6 

90-15.6E-6 

120-1.59 

120-1.21 

120-  9.5E-6 

120-7.9E-6 


I  'Ml 

f 


r* 


IK* 


t 


iii.- 

wm 
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MD  2005  — Record  49  Pistol  280  m 


160-0.45 

160-0.49 

160-1.2E-6 

160-3.7E-6 

100-1.06 

100-0.67 

100-6.0E-6 

100-5.2E-6 

40-3.37 

40-1.87 

40-  3.9E-6 

40-  8.6E-6 

60-11.5 

60-2.03 

60-17. 9E-6 

60-14. OE-6 

90-8.50 

90-1.91 

90-13.4E-6 

90-15. IE-6 

120-3.08 

120-1.29 

120-  9. OE-6 

120-  8.7E-6 


P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 


-0.5  0  0.5  1  1.5  2 


Time,  sec 
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MD  2005  —  Record  50  Pistol  280  m 


160-0.58 

160-0.30 

160-  1.3E-6 

160-  1.9E-6 

100-0.84 

100-0.50 

100-4.6E-6 

100-4.7E-6 

40-3.34 

40-1.74 

40-  4.5E-6 

40-  8.6E-6 

60-13.1 

60-2.02 

60-18.2E-6 

60-12.9E-6 

90-7.77 

90-1.66 

90-11.7E-6 

90-13.5E-6 

120-1.95 

120-1.10 

120-8.1E-6 

120-  7.7E-6 


P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 


-0.5  0  0.5  1  1.5  2 


Time,  sec 
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MD  2005  — Record  51  Pistol  280  m 


160-0.42 

160-0.34 

160-  1.8E-6 

160-2.5E-6 

100-0.54 

100-0.53 

100-5.6E-6 

100-4.6E-6 

40-2.33 

40-1.45 

40-  5.5E-6 

40-  8.2E-6 

60-8.49 

60-1.82 

60-17. 2E-6 

60-11.7E-6 

90-6.03 

90-1.61 

90-11.3E-6 

90-12.7E-6 

120-2.00 

120-1.06 

120-8.0E-6 

120-8.5E-6 


P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 


-0.5  0  0.5  1  1.5  2 


Time,  sec 
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MD  2005  —  Record  52  Pistol  280  m 


160-0.51 
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160-  3.2E-6 

100-1.13 

100-0.75 

100-6.4E-6 

100- 5. IE-6 

40-2.93 

40-1.65 

40-  4.1  E-6 

40-  8.4E-6 

60-10.7 

60-1.85 

60-17. IE-6 

60-12. 9E-6 

90-6.56 

90-1.70 

90-11.8E-6 

90-14.4E-6 

120-1.88 

120-1.24 

120-8.8E-6 

120- 9. IE-6 


rtf* 


p 

v-' 

R 


4~ 
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"ipp’rr 

_4_ 


l* 


IK 
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p- 
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MD  2005  —  Record  53  Pistol  310  m 


190-0.38 

190-0.30 

190-5.0E-6 

190-5.4E-6 

130-0.38 

130-0.43 

130-  3.8E-6 

130-  3.9E-6 

70-1.05 

70-0.62 

70-2.0E-6 

70-3.3E-6 

30-22.6 

30-4.91 

30-47. OE-6 

30-21 .7E-6 

60-10.4 

60-3.01 

60-23.4E-6 

60-22.4E-6 

90-2.74 

90-1.63 

90-13.7E-6 

90-  9.6E-6 


P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 


-0.5  0  0.5  1  1.5  2 


Time,  sec 
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MD  2005  — Record  54  Pistol  310  m 


190-0.37 

190-0.31 

190-1.9E-6 

190-2.4E-6 

130-0.57 

130-0.42 

130-4.3E-6 

130-3.2E-6 

70-1.05 

70-0.80 

70-1.8E-6 

70-  4.1  E-6 

30-36.0 

30-8.33 

30-53.1  E-6 

30-26.1  E-6 

60-20.4 

60-4.70 

60-27. OE-6 

60-27.1  E-6 

90-5.42 

90-2.00 

90-15.9E-6 

90-13. OE-6 


P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 


-0.5  0  0.5  1  1.5  2 


Time,  sec 
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MD2005  —  Record  55  Pistol  310  m 


190-0.34 

190-0.38 

190-  1.4E-6 

190-  2.3E-6 

130-0.52 

130-0.59 

130-  4.7E-6 

130- 4. IE-6 

70-1.33 

70-0.95 

70-  1.7E-6 

70-  4.3E-6 

30-43.1 

30-9.18 

30-58.1  E-6 

30-30. 9E-6 
60-18.5 

60-4.07 

60-30. 7E-6 

60-30.3E-6 

90-4.17 

90-2.03 

90-16. IE-6 

90-14.3E-6 


P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 


-0.5  0  0.5  1  1.5  2 


Time,  sec 
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MD  2005  —  Record  56  Pistol  310  m 


190-0.35 

190-0.32 

190-  1.4E-6 

190-3.2E-6 

130-0.57 

130-0.49 

130-4.4E-6 

130-3.8E-6 

70-1.36 

70-0.83 

70-1.9E-6 

70-  4.2E-6 

30-27.5 

30-6.26 

30-46.7E-6 

30-24.5E-6 

60-11.9 

60-3.00 

60-24. 5E-6 

60-23. 2E-6 

90-2.13 

90-1.71 

90-13.8E-6 

90-11.6E-6 


P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 


-0.5  0  0.5  1  1.5  2 


Time,  sec 
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MD  2005  — Record  57  Pistol  310  m 


190-0.39 

190-0.36 

190-  1.6E-6 

190-2.9E-6 

130-0.45 

130-0.44 

130-  3.6E-6 

130- 4. IE-6 

70-1.15 

70-0.84 

70-  1.8E-6 

70-4.9E-6 

30-41.7 

30-7.88 

30-70. OE-6 

30-37. 6E-6 

60-17.7 

60-5.00 

60-37.3E-6 

60-33.8E-6 

90-3.49 

90-2.52 

90-20. 3E-6 

90-15.7E-6 


P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 


-0.5  0  0.5  1  1.5  2 


Time,  sec 
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MD  2005  — Record  62  Pistol  430  m 


310-0.48 

310-0.32 

310-  1.5E-6 

310-2.3E-6 

250-1.06 

250-0.51 

250-  3.6E-6 

250-  3.9E-6 

190-1.19 

190-0.68 

190-  1.9E-6 

190-2.9E-6 

90-2.24 

90-1.49 

90-14.5E-6 

90-7.2E-6 

60-7.14 

60-2.46 

60-27.0E-6 

60-9.2E-6 

30-39.8 

30-8.70 

30-54. OE-6 

30-31 .4E-6 


P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 


-0.5  0  0.5  1  1.5  2 


Time,  sec 
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MD  2005  — Record  64  Pistol  280  m 


160-0.86 

160-0.71 

160-1.5E-6 

160-6.2E-6 

100-4.08 

100-1.34 

1 00-1 3.5E-6 

100-1 1.8E-6 

40-22.4 

40-4.99 

40-14. IE-6 

40-22.3E-6 

60-0.43 

60-0.45 

60-4.9E-6 

60-4.0E-6 

90-0.35 

90-0.42 

90-4.8E-6 

90-3.9E-6 

120-0.22 

120-0.32 

120-2.5E-6 

120-2.3E-6 


P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 

P-340-1 

P-340-0 

V-340 

R-340 

P-370-1 

P-370-0 

V-370 

R-370 

P-400-1 

P-400-0 

V-400 

R-400 


-0.5  0  0.5  1  1.5  2 


Time,  sec 
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MD  2005  —  Record  107  Cap  P+2+2 


PP-10-0 

PP-22-1 

PP-32-1 

PP-32-0 

PP-62-1 

PP-62-0 

VP-62 

RP-62 

PP-92-0 

PP-92-1 

VP-92 

RP-92 

PB-10-0 

PB-22-1 

PB-32-1 

PB-32-0 

PB-62-1 

PB-62-0 

VB-62 

RB-62 

PB-92-1 

PB-92-0 

VB-92 

RB-92 

P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 
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MD  2005  —  Record  109  Cap  BP+2+2 


PP-10-0 

PP-22-1 

PP-32-1 

PP-32-0 

PP-62-1 

PP-62-0 

VP-62 

RP-62 

PP-92-0 

PP-92-1 

VP-92 

RP-92 

PB-10-0 

PB-22-1 

PB-32-1 

PB-32-0 

PB-62-1 

PB-62-0 

VB-62 

RB-62 

PB-92-1 

PB-92-0 

VB-92 

RB-92 

P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 
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MD  2005  —  Record  110  Cap  BP+2+2 


PP-10-0 

PP-22-1 

PP-32-1 

PP-32-0 

PP-62-1 

PP-62-0 

VP-62 

RP-62 

PP-92-0 

PP-92-1 

VP-92 

RP-92 

PB-10-0 

PB-22-1 

PB-32-1 

PB-32-0 

PB-62-1 

PB-62-0 

VB-62 

RB-62 

PB-92-1 

PB-92-0 

VB-92 

RB-92 

P-120-1 

P-120-0 

V-120 

R-120 

P-180-1 

P-180-0 

V-180 

R-180 

P-240-1 

P-240-0 

V-240 

R-240 
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MD  2005 --Record  111  Cap  BP+2+2 


Time,  sec 
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Appendix  A:  MATLAB  Files  for  Blast  Absorber  Test 

The  binary  data  files  from  the  NZ  seismograph  were  analyzed  and  plotted  using 
the  Matlab  files  listed  in  this  section.  These  files  include: 


doplotMD05.m  -  this  file  sets  the  record  numbers  to  process,  assigns  the  directory 
where  the  data  files  are  stored  (something  that  will  have  to  be  changed  by  the 
user),  reads  in  the  SEG2  format  binary  data,  converts  voltages  to  physical  units, 
and  plots  the  data.  This  program  calls  the  following  subroutines: 

readnzATC.m  -  read  the  SEG-2  binary  data  using  the  field  record  numbers 

words.m  -  parse  ascii  header  lines  into  individual  words 

plotNZ48FourCol.m  -  plot  the  C-4  data  in  four  panels  OR 

plotNZ48.m  -  to  plot  the  pistol  or  blasting  cap  data  on  one  panel 

MD051abel.m  -  to  get  sensor  geometry,  sensitivities,  etc 


The  conversion  to  physical  units  is  done  in  the  plotxx.m  subroutine  using  the  in¬ 
formation  from  the  MD051abel.m  subroutine. 

doplotMD05.n1 

%  Fujitsu  C:\aaaNZData\PlotATC05\DataReportFiles\doplotMD05.m 

%  Reads  NZ  binary  records,  makes  3  panel  plots  of  data 
%  Calls  functions 

%  readnzATC.m  -  to  read  the  SEG-2  binary  data  using  the 
%  field  record  numbers 

%  words.m  -  to  parse  ascii  header  lines 
%  plotnzmout.m  -  to  plot  the  data 

%  MD051abel.m  -  to  get  sensor  geometry,  sensitivities,  etc 

%  NZ  binary  record  numbers; 

recs  =  [106:111];  %Blasting  caps  NOT  48  channel 

recs  =  [48:67];  %A11  pistol  shots 

recs  =  [37:40  42:46  69:92  ];  %A11  C4  shots 

%  Loop  to  plot  the  files 
for  i  =  1 : length ( recs ) ; 
recnum  =  recs (i)  ; 
shotnum  =  recs (i) ; 

%  Construct  binary  data  file  name 
fname  =  sprintf ( . . . 
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1 C : \\DDrive\\aaaFY05BlastAbsorberTest\\BlastAbsorberCRRELBinaryDa 
ta\\%s . dat ' , . . . 

num2str (recnum) ) ; 


%  Read  in  data  from  NZ  binary  file 

[x, npts , deltat , nchan, delayms , descalingf act , . . . 

stackcount ,  acqdate, acqtime]  =  readnzATC (fname, shotnum) ; 


%  Plot  the  binary  data 

%  plotNZ48;  %  Use  for  pistol  shots  and  blasting  caps 

plotNZ48FourCol;  %  For  C4  shots 
end  %loop  over  files 


return 


readnzATC.m  -  read  the  SEG-2  binary  data  using  the  field  record  numbers 

This  function  reads  a  binary  seismograph  file  in  SEG-2  format,  and  returns  the 
data  and  header  variables  to  the  calling  program.  It  is  called  by  doplotMD05.m 
and  does  not  call  any  other  subroutines. 


function  [ scan ,  samplesPerScan, samplinglnterval , nbOf Traces ,  ... 

delaytime, descalingf act , stackcount ,  acqdate, acqtime]  =  ... 
readnzATC  ( filename, shotnum) 

%This  version  uses  shotnum  to  select  correct  STACK  index  in 
header 

%  Lines  158  below 

%  This  fen  reads  a  binary  seismograph  file  in  SEG-2  format,  and 
returns  the 

%  data  in  "scan"  variable  above. 

%  This  code  written  by  Don  Albert,  CRREL,  and  based  on  code 
"SEG2L0AD"  written 

%  by  Piece  PY  to  read  a  radar  file  in  modified  SEG-2  format. 
Some  of  Piece  1 s 

%  original  code  is  still  here  in  this  file. 

%  Calls  fen  words. m 

%  To  read  one  NZ  file  directly: 

%fname  = 

1 10077 . dat '  ; [x, npts , deltat, nchan, npts 2 , delayms , descalingf act , acqd 
ate, acqtime]  =  readnz (fname) ; 

%  READNZ  -  For  Geometries  NZ  seismograph 

%  READNZ  Modified  by  D  Albert  to  read  standard  SEG-2  format 
%  SEG2LOAD  Read  a  SEG-2  (standard  SEG-2  format  of  the  Society  of 
Exploration  Geophysicist) 

%  file  from  disk. 

o 

o 

[scan, samplesPerScan, samplinglnterval, shaf tlnterval , timerFrequenc 

y]  = 

%  =seg21oad  ('filename')  reads 
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%  the  file  'filename'  and  returns  the  image  scan  [m,n] 

containing 

%  n  A-Scan  of  m  samples. 

%  If  no  extension  is  given  for  the  filename,  the  exten¬ 

sion  ' . sg2 '  is  assumed. 

%  samplesPerScan  contains  the  number  of  samples  per  A- 


scan 

%  shaf tlnterval  contains  the  distance  between  shaft 

encoder  triggers  in  meter 

%  samplinglnterval  contains  the  time  between  2  samples 

in  pico-seconds 

%  timerFrequency  contains  the  frequency  of  A-scan  sam¬ 


pling  in  Hz 

o 

o 


%  Piece  PY  24/07/1996 

%  LAMI  -  DeTeC  Demining  Technology  Center 

%  Swiss  Federal  Institute  of  Technology  (EPFL)  - 

Lausanne,  Switzerland 


%  check  argument  and  filename 
if  (nargin==0) 

error  ( ' readseg2  requires  at  least  a  filename  as  an  argument 

!  ')  ; 

end; 

if  (isstr  ( filename) ~=1 ) 

error  ('Argument  is  not  a  filename  !'); 
end; 

if  (isempty  (findstr  ( filename, '.')) ==1 ) 
filename  =[ filename, '. sg2 '] ; 
end; 

f id=f open  (filename, ' rb ' , ' ieee-le ' ) ; 
if  (fid  ==-l) 

error  (['Error  opening  ', filename,'  for  input  !']); 
end; 


%  check  for  SEG-2  file  type 

%  first  2  bytes  equal  '3a55h'  (14933)  for  PC/Windows 

f ileType=f read  ( fid, 1 ,' short ') ; 
if  (fileType  ~=  14933) 
fclose  (fid) ; 

error  ('Not  a  SEG-2  file  !'); 
end; 


%  Open  a  log  file  to  keep  track  of  files  that  were  read 
fid2  =  f open  (  ' readnzlog . m ' ,  ' a ' )  ; 
fprintf ( f id2 ,  '  \n'); 

samplesPerScan  =0; 
shaf tlnterval  =0; 
samplinglnterval=0 ; 
timerFrequency  =0; 

%  read  File  Descriptor  Block 

revNumber  =  fread  ( fid, 1 ,' short ') ; 

sizeOf TracePointer  =  fread  ( f id,  1 ,  ' ushort ' ) ; 
nbOfTraces  =  fread  ( fid, 1 ,' ushort ') ; 


o\°  o\° 
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sizeOf ST 

=  fread 

(fid, 1 , ' uchar 1 ) ; 

f irstST 

=  fread 

(fid, 1, 'char' ) ; 

secondST 

=  fread 

(fid, 1, 'char' ) ; 

sizeOf LT 

=  fread 

(fid, 1 , ' uchar ' ) ; 

firstLT 

=  fread 

(fid, 1, 'char' ) ; 

secondLT 

=  fread 

(fid, 1, 'char ' ) ; 

reserved 

=  fread 

( f id, 18 , ' uchar ' ) ; 

tracePointers 

=  fread 

(fid, nbOf Traces , ' ulong ' ) ; 

%  read  free  strings 

f seek  ( f id, 32  +  sizeOf TracePo inter ,  '  bof ' ) ; 
offset  =  fread  ( f id,  1 ,  ' ushort ' ) ; 

%  File  descriptor  block 
while  (offset  >  0) 

freeString  =  setstr  (fread  ( fid, of f set-2 ,  1  char ' ))  1 ;  % 

%Decode  file  descriptor  block 

if  (findstr  ( freeString, ' ACQUIS ITION_DATE 1 )  >  0) 

acqdate  =  (freeString  . . . 

(length  (' ACQUIS IT ION_DATE  ') : length  (freeString))); 

end 

if  (findstr  ( freeString,  1  ACQUIS ITION_TIME ' )  >  0) 

acqtime  =  (freeString  . . . 

(length  (' ACQUIS IT ION_TIME  ') : length  (freeString))); 

end 

offset  =  fread  ( fid, 1 ,  ' ushort ') ; 
end; 

%  acqdate,  acqtime,  offset 
%  read  traces 

find  number  of  samples  per  trace 
First  trace  descriptor  block 
fseek  ( f  id,  tracePointers  (l),'bof'); 
traceld  =  fread  (fid, 1, 'ushort ') ; 

sizeOfBlock  =  fread  ( fid, 1 ,' ushort ')  ; 
sizeOfData  =  fread  ( f id, 1 , ' ulong 1 ) ; 
nbOfSamples  =  fread  ( fid, 1 ,' ulong ') ; 
samplesPerScan  =  nbOfSamples; 

%  traceld,  sizeOfBlock,  sizeOfData, nbOfSamples, samplesPerScan 

dataCode  =  fread  ( f id,  1 ,  ' uchar 1 ) ; 
reserved  =  fread  ( fid, 19,  * uchar 1 ) ; 
offset  =  fread  (fid, 1, 'ushort ') ; 
while  (offset  >  0) 

freeString  =  setstr  (fread  (fid, of f set-2 ,' char '))' ; 
if  (findstr  ( freeString, ' SAMPLE_INTERVAL ' )  >  0) 

samplinglnterval  =  str2num  (freeString  . . . 

(length  ( ' SAMPLE_INTERVAL  ' ) : length  (freeString) ) ) ; 

end 

if  (findstr  ( freeString, ' DELAY ' )  >  0) 

delaytime  =  1000*str2num  (freeString  . . . 

(length  ('DELAY  '): length  (freeString))); 

end 

if  (findstr  ( freeString, ' DESCALING_FACTOR ' )  >  0) 
descalingf actor  =  str2num  (freeString  . . . 
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(length  ( 1 DESCALING_FACTOR  '): length  (freeString))); 

end 

offset  =  fread  ( f id, 1 , ' ushort 1 ) ; 
end; 

%  samplinglnterval ,  delaytime,  descalingf actor 

%Write  file  descriptor  to  log  file 

fprintf ( f id2 , ' %s  %s  %s  %g  %g  %g  %g 
\n  ’  , filename, acqdate, acqtime, nbOf Samples , . . . 

samplinglnterval, delaytime, descalingf actor ) ; 

scan  =  zeros  (nbOf Samples , nbOf Traces )  ; 
for  i=l : nbOf Traces , 

fseek  ( f id, tracePointers  (i),'bof'); 
traceld  =  fread  (fid, 1, 'ushort 1 ) ; 
sizeOfBlock  =  fread  ( fid, 1 ,' ushort ')  ; 
sizeOfData  =  fread  ( f id, 1 , ' ulong 1 ) ; 
nbOfSamples  =  fread  ( fid, 1 ,' ulong '  ) ; 
dataCode  =  fread  ( f id, 1 , ' uchar 1 ) ; 
reserved  =  fread  (fid, 19, 1 uchar '  ) ; 

%  Trace  descriptor  blocks 
%Reads  all  at  once 

freeString  =  setstr  (fread  ( fid, sizeOfBlock-32 ,' char '))’ ; 
w  =  words ( freeString) ; 

%  Remove  percent  below  to  print  trace  header  to  screen  for 
%  adjustment 
%if  ( i==l )  ;w,end 

descalingf act ( i )  =  str2num (w (7, 1 : 14) ) ; 
if (shotnum  <  10  |  shotnum  >  99) 

stackcount ( i )  =  str2num (w (27 , 1 : 2 ) ) ;  %files  1  -  9.dat  and 
>=100 . dat 
else 

stackcount ( i )  =  str2num (w (28 , 1 : 2 ) ) ;  %files  >=10. dat  and 

<=99 . dat 
end 

%  w  =  words ( freeString) ; 

%MOUT  data 

%descalingf act ( i )  =  str2num (w (4, 1 : 14) ) ; 

%stackcount (i)  =  str2num (w ( 18 , 1 : 2 ) ) ; 

%  descalingf act ( i )  =  str2num (w (7 , 1 : 14 ) ) ; 

%stackcount ( i )  =  str2num (w (2 6, 1 : 2 ) ) ; 

%  stackcount ( i )  =  str2num (w (27 , 1 : 2 ) ) ; 

%  read  the  data  here 

scan  ( 1 : nbOfSamples , i )  =  fread  ( fid, nbOfSamples ,' float32 ') ; 

end; 

fclose  (fid) ;  %close  data  file 
return 


%  CHANNEL_NUMBER 
%  48  DELAY 

%  -0.050  !  DESCALING_FACTOR 
%  2 . 697400E-003  FIXED  GAIN 
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%  24 
%  DB 

%  LINE_ID 
%  0  LOW_CUT_FILTER 
%  0 

%  0  NOTCH_FREQUENCY 
%  0  RAW_RECORD 
%  RECEIVER_LOCATION 
%  570.00  SAMPLE_INTERVAL 
%  0.000125  SKEW 
%  -0.0000625  SOURCE_LOCATION 
%  0.00 
%  STACK 
%  1  NOTE 
%  DISPLAY_SCALE 
%  67 

%%SKEW  =  ??? 


%  MN02 : 

%  ALIAS_FILTER 
%3333 . 33 

%0  AMPLITUDE_RECOVERY 
%NONE  CHANNEL_NUMBER 
%1  DELAY 

%-0 . 050  !  DESCALING_FACTOR 

%2 . 697400E-003  -  DIGITAL_HIGH_CUT_FILTER 

%0 

%0  D I G I TAL_LOW_CUT_F I LTER 

%0 

%0  FIXED_GAIN 

%24 

%DB 

%  LINE_ID 

%0  LOW_CUT_FILTER 
%0 

%0  NOTCH_FREQUENCY 
%0  RAW_RECORD 

%D : \4000\35 . dat  RECEIVER_LOCATION 
%510 . 00  SAMPLE_INTERVAL 
%0. 000125  SHOT_SEQUENCE_NUMBER 
%35  SKEW 

%-0. 0000625  SOURCE_LOCATION 

%0.00 

%  STACK 

%1  NOTE 

%DISPLAY_SCALE 

%63 

words.m  -  parse  ascii  header  lines  into  individual  words 


This  function  separates  a  long  string  of  text  into  separate  words.  It  does  not  call 
any  other  subroutines.  This  function  is  from  the  Matlab  manual. 

function  all_words  =  words ( input_string) 

%  filename:  words.m 
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%  Separates  words  in  a  long  string 
%  Individual  words  are  then  in  all_words ( 9 , : ) 
%  From  the  book.  Using  Matlab,  p.  11-12 

remainder  =  input_string; 
all_words  =  ' ' ; 

while  (any (remainder ) ) 

[ chopped, remainder ]  =  strtok ( remainder ) ; 
all_words  =  strvcat (all_words , chopped) ; 

end 

return 


plotNZ48FourCol.m  -  plot  the  C-4  data  in  four  panels 

This  code  plots  48  channels  of  data  in  four  separate  panels  on  a  sheet  of  landscape 
oriented  paper.  Because  of  intermittent  interference  with  the  blasting  electrical 
pulse,  the  program  mutes  or  zeros  the  first  20  ms  of  data  on  channels  1  -  24.  In 
addition,  the  first  226  ms  of  data  for  channel  15  are  muted  since  that  channel  suf¬ 
fered  from  additional  interference. 

This  code  calls  MD051abel.m  to  get  the  experimental  information  including 
source  type  and  location,  receiver  location,  receiver  sensitivities,  etc. 

%  PlotNZ4 8FourCol . m  -  48  channels  in  4  columns 
%  Includes  TB  muting  for  channels  1-24 

%  Need  to  set  file  name  and  read  in  data  first 
%  using  readnz  sbr 
%  Calls  sbrs 

%  MD051abel.m  -  supplies  the  details  of  the  expt 
%  including  sensitivity,  coords  of  each  channel  etc. 

%  Don  Albert  Donald.G.Albert@erdc.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 

iprint=l;  %  to  print  and  delete  plot 
iprint  =  0;  %skip  printing 

nchan  =  48; 

srate  =  1/  deltat;  %number  of  samples  per  second 
t  =  ( 1 : npts ) *deltat  +  delayms/1000 ;  % 
tend  =  npts; 

t start 2=1 ; tend2=tend; t2=t ; 
tstartpanell  =  0.0; 
tstartpanell  =  delayms/1000;  % 
tendpanell  =  tstartpanell  +  t (tend) ; 

%  Default  time  window  (entire  record) 

tt2  =  0.0;  tt3  =  t(tend); 

tt2  =  tstartpanell;  tt3  =  t (tend) ;  % 


ERDC/CRREL  TR-07-3 


74 


tstart2  =  f loor (tt2/deltat ) ; tend2  =  round (tt3/deltat) ; 
if  (tstart2  <=  1);  tstart2  =  1;  end 
tstartpanell  =  tt2 ; tendpanell  =  tt3; 

%apply  TB  mute  to  the  close-in  pressure  sensors  and  VB-62  (Chan 
18) 

imutechannels  =  [3  4  7:12  15  16  18:24] ; 

%  Apply  TB  mute  to  the  close-in  channels  (geos  infrequently  need 
it) 

imutechannels  =  [1:24] ; 
tmute  =  0.02; 

indexmute  =  round (tmute/deltat) -delayms/ ( 1000*deltat) ; 
x ( 1 : indexmute, imutechannels )  =  0; 

%  PB-92-0  needs  add'l  muting 

tmute  =  0.226; 

indexmute  =  round (tmute/deltat ) -delayms/ ( 1000*deltat) ; 
x ( 1 : indexmute, 15 )  =  0; 

xx  =  [  ] ; 

for  i  =  1 : nchan; xmax (i) =max (abs (x ( : , i) ) ) ;  end 

for  i  =  1 : nchan; xmax (i) =max (abs (x (tstart2 : tend2 , i) )) ;  end 

if (xmax (i)  ==  0);  xmax(i)  =  1;  end 

%  Get  geometry,  etc  for  this  experiment 

[sensor, sensortype, rtype, sensorx, sensory, sensorz,  .  .  . 

srcloc,  sourcetype, stype, srcx, srcy, srcz, sourcesize,  .  .  . 
engf act , dist , sensor2 , source2 ]  =  MD051abel 
(x, shotnum, nchan) ; 

if ( shotnum>=l 12  &  shotnum  <=117 ) ;  nchan  =  36;  end 
if (shotnum>=118  );  nchan  =  24;  end 

%  Conversion  factor  from  voltage  to  physical  units  (Pa  or  m/s) 
%xgain  =  descalingf act . / (stackcount . *engf act ) ;  %recorder  units  to 
eng  units 
for  i  =  1: nchan; 

if ( stackcount ( i ) ==0 ) ;  stackcount (i) =  1 ;  end 
%recorder  units  to  eng  units 

xgain(i)  =  descalingf act ( i )./( stackcount ( i ). *engf act ( i ))  ; 
xx  (  : ,  i)  =  x(:,i)  *  xgain(i); 

%xxmax(i)  =  max (abs (xx) ) ;  %abs  max  for  plotting  -  entire  trace 
xxmax(i)  =  max (abs (xx (tstart2 : tend2 , i) ) ) ;  %abs  max  for  plot¬ 
ting 
end 

[aa  bb]  =  max (xx) ;  %set  plot  windows  because  manual  TB 
[cc  dd]  =  min (xx) ;  %for  output  file 

tmaxl  =  delayms/1000  +  (bb(l)  ) /srate;  %P-l-0  times  in  seconds 

tmax21  =  delayms/1000  +  (bb (21 )) /srate;  %P-100-0 

tstartl=  floor (tmaxl  *10) /10; 

tstart21=  f loor (tmax21*10 ) /10 ; 

textmax=max (max ( t ) ) +max (max ( t ) ) / 3  0 ; 

textmin=min ( t ) -max (max ( t ) ) / 3  0 ; 

nametext  =  [*NZ  Data  -  1  panel'];  % - FIGURE  -  4 

Panel 

fl  =  f igure (' Name ', nametext ) ; 
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set  (gcf ,  ' Units ' ,  1  inches ’ ,  ' PaperUnits 1 ,  1  inches ' ) ; 
iplot  =  [37:48  25:36  13:24  1:12  ]; 

%  Switch  channels  close  to  source  for  bare  ground  shots: 
if  (srcloc (1 ) == 1 B 1 ) ; 

iiplot  =  [37:48  25:36  1:12  13:24  ]; 

end 

%  Use  same  channel  order  for  all  shots 
iiplot  =  [37:48  25:36  1:12  13:24  ]; 


ooooooooooooooooooooooooool 


subplot (1,4,1) ; 
tstartpanell=0 . 0 ;  tendpanell=0 . 4 ; 
%  First  12  channels 
iplot=iiplot (37 : 48)  ; 
nplot  =  length (iplot ) ; 


o  o  o  o 
o  o  o  o 


%  Set  axis  params  here;  depends  on  nplot  -  t2 
yshift  =  (l:nplot)*2  -  2; 

%ax  =  [t2 (1)  max(t2)  -1  2*nplot] ; 

ax  =  [ tstartpanell-0 . 05  tendpanell  -1  2*nplot] ; 

%ax  =  [13-1  2*nplot] ; 
axis  (ax) ; 

axl  =  gca;  set  (axl ,  1 YTick ’ ,  [ ] ) ;  set (axl ,  ' Box  1 ,  ' on ' ) ; 
xlabel ( ' Time,  sec1) 


%  plot  data  and  label 

hold  on; 

for  i=l : nplot 

plot(t,  (xx (1 : length (t) , iplot (i) ) /xxmax (iplot (i) ) ) 
+yshift (i) , 1 k! ) ; 

leftlabel  =  sprintf (' %g '  , round (dist (iplot (i) )))  ; 
distlabel  =  sprintf ( 1 %g ', round (dist (iplot (i) ))) ; 
strl  =  sensor ( iplot ( i ),:)  ; 

str2  =  sprintf ( ' %7 . 2e 1 , (xxmax ( iplot ( i ) ) ) ) ; 
if  (xxmax ( iplot ( i )) )  <  0.0001 

str2=sprintf ( 1 %4 . lfE-6 ' , (lE6^xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 

str2=sprintf ('%4.2fE-3', ( IE 3* xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.01 

str2=sprintf ( 1 %4 . 2  f E-3 1 ,  ( lE3^xxmax ( iplot ( i ) ) ) )  ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf ( ' %4 . lfE-3 1 , (lE3^xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot  ( i )) )  <  10.0 

str2=sprintf ( 1 %4 . 2  f ' ,  (xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot  ( i )) )  <  100.0 

str2=sprintf ( ' %4 . If ' , (xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot  ( i )) )  >  1000.0 

str2=sprintf ( ' %4 . If k ' , (  xxmax (iplot (i))/1000  )); 
else 

str2=sprintf ( ' %4 . If ' ,  (xxmax (iplot (i) ) )  )  ; 
end  %text  loop 
str3  =  char (strl); 

str4  =  [str3  str2];  %Concat  on  one  line 

str4  =  [str3  ] ;  %Concat  on  one  line 

str5  =  cellstr (str4 ) ; 

lef tlabel02  =  [leftlabel  str2]  ; 

str4  =  char(str2  ,  str3) ; 
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str5  =  cellstr ( str4 ) ; 

text (tstartpanell-0 .05-0.02,  [y shift (i) ] , str5,  ' Font Size '  ,  8,  .  .  . 

' HorizontalAlignment 1 ,  1  right ' ) ; 

text ( tendpanell-0 .02, [y shift (i) +0 . 3] , dist label, ' Font Size ' , 8,  .  .  . 
' HorizontalAlignment ' ,  1  right ' ) ; 

end 


Qnhnl  n+-  (  1  A  ?)  •  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  9-  ' 

oUJjpiUL  ^1, 4,  Z,/  ,  ooooooooooooooooo 

tstartpanell=0 . 0 ;  tendpanell=0 . 4 ; 

%  First  12  channels 


o  o  o 
o  o  o 


iplot=iiplot (25 : 36)  ; 
nplot  =  length (iplot )  ; 

%  Set  axis  params  here;  depends  on  nplot  -  t2 
yshift  =  (l:nplot)*2  -  2; 

ax  =  [ tstartpanell-0 . 05  tendpanell  -1  2*nplot] ; 
axis (ax) ; 

axl  =  gca;  set (axl ,  1 YTick 1 ,  [ ] ) ;  set (axl ,  ' Box  1 ,  ' on ' ) ; 
xlabel ( ' Time,  sec  1 ) 

%  plot  data  and  label 
hold  on; 


for  i=l : nplot 

plot (t, (xx (1 : length (t) , iplot (i) ) /xxmax (iplot (i) ) ) 
+yshif t ( i ) , 1 k 1 ) ; 

leftlabel  =  sprintf (’ %g 1 , round (dist (iplot (i) ))) ; 
distlabel  =  sprintf (' %g ', round (dist (iplot (i) ))) ; 
strl  =  sensor ( iplot ( i ),:) ; 

str2  =  sprintf ( 1 %7 . 2e ' , (xxmax ( iplot ( i ) ) ) ) ; 
if  (xxmax ( iplot ( i )) )  <  0.0001 

str2=sprintf ( 1  %  4 . lfE-6 ' ,  (lE6^xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 

str2=sprintf ( 1 %4 . 2  fE-3 1 ,  ( IE 3* xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.01 

str2=sprintf ( *  %4 . 2 f E-3 1 ,  ( IE 3^ xxmax ( iplot ( i ) ) ) )  ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf ( ' %4 . lfE-3 ' , (lE3*xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  10.0 

str2=sprintf ( ' %4 . 2  f ' ,  (xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i ))  )  <  100.0 

str2=sprintf ( ' %4 . If ' , (xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  >  1000.0 

str2=sprintf ( ' %4 . If k 1 , (  xxmax (iplot (i))/1000  )); 
else 


str2=sprintf ( ’ %4 . If ' , (xxmax (iplot (i) ) ) ) ; 
end  %text  loop 
str3  =  char (strl); 
str4  =  char(str2  ,  str3)  ; 
str5  =  cellstr ( str4  )  ; 

text  (tstartpanell-0 .05-0.02,  [yshift (i) ] , str5,  ' Font Size  1 , 8,  .  .  . 

' HorizontalAlignment ' , ' right ' ) ; 

text (tendpanell-0 .02,  [yshift (i) +0 . 3] , distlabel,  1  Font Size ' , 8,  .  .  . 
' HorizontalAlignment ' ,  ' right 1 ) ; 


subplot (1,4,3) 
tstartpanell=0 


2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-9-Prin^l  ^ 

ooooooooooooooooooooooooo  o  r  cLllt^  _L  O  oooooooo 


3;  tendpanell=0 . 9 ; 
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%  120-240  m 
iplot=iiplot(13:24) ; 
nplot  =  length (iplot) ; 

%  Set  axis  params  here;  depends  on  nplot  -  t2 
yshift  =  (l:nplot)*2  -  2; 

ax  =  [ tstartpanell-0 . 05  tendpanell  -1  2*nplot] ; 
axis (ax) ; 

axl  =  gca;  set  (axl ,  1 YTick ’ ,  [ ] ) ;  set (axl ,  ' Box ' ,  ' on ' ) ; 
xlabel ( ' Time ,  sec') 

%  plot  data  and  label 

hold  on; 

for  i=l : nplot 

plot  (t,  (xx (1 : length (t) , iplot (i) ) /xxmax (iplot (i) ) ) 

+yshif t ( i ) , ' k 1 ) ; 

leftlabel  =  sprintf (' %g '  , round (dist (iplot (i) )))  ; 
distlabel  =  sprintf (' %g ', round (dist (iplot (i) ))) ; 
strl  =  sensor ( iplot ( i ),:)  ; 

str2  =  sprintf ( 1 %7 . 2e 1 ,  (xxmax ( iplot ( i ) ) ) ) ; 
if  (xxmax ( iplot  ( i )) )  <  0.0001 

str2=sprintf ( 1  % 4 . lfE-6 ' ,  (lE6*xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 

str2=sprintf ( 1 %4 . 2 fE-3 1 ,  ( IE 3* xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.01 

str2=sprintf ( ' %4 . 2  fE-3 ' ,  ( IE 3* xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf ( 1 %4 . lfE-3 ' ,  (lE3*xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  10.0 

str2=sprintf ( 1 %4 . 2  f ' ,  (xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  100.0 

str2=sprintf ( '%4 . If' ,  (xxmax (iplot (i) )  )  )  ; 
elseif  (xxmax ( iplot ( i ))  )  >  1000.0 

str2=sprintf ( 1 %4 . If k 1 , (  xxmax (iplot (i))/1000  )); 
else 

str2=sprintf ( 1 %4 . If ' ,  (xxmax (iplot (i) ) )  )  ; 
end  %text  loop 
str3  =  char (strl); 
str4  =  char(str2  ,  str3)  ; 
str5  =  cellstr ( str4 ) ; 

text (tstartpanell-0 .05-0.02, [yshift (i) ] , str5, ' Font Size ' , 8, . . . 

' HorizontalAlignment ' ,  1  right ' ) ; 

text ( tendpanell -0 .02, [yshift (i) +0 . 3] , distlabel, ' Fonts ize ' , 8, . . . 
1 HorizontalAlignment 1 ,  ' right 1 ) ; 

trllUL  oooooooooooooooooooooooooooooooooooooooooooooooooooooo 


subplot (1,4,4)  ; 
tstartpanell=0 . 9 
%  340-400  m 


oooooooooooooooooooooooooo 


;  tendpanell=l . 4 ; 


Panel 


4 


oooooooo 


iplot=i iplot (1:12)  ; 
nplot  =  length (iplot) ; 

%  Set  axis  params  here;  depends  on  nplot  -  t2 
yshift  =  (l:nplot)*2  -  2; 

ax  =  [ tstartpanell-0 . 05  tendpanell  -1  2*nplot] ; 
axis (ax) ; 

axl  =  gca;  set (axl ,  f YTick 1 ,[]) ;  set (axl ,  ’ Box  1 ,  ' on  1 ) ; 
xlabel (' Time,  sec') 

%  plot  data  and  label 
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hold  on; 
for  i=l : nplot 

plot  (t,  (xx (1 : length (t) , ip lot (i) ) /xxmax (iplot (i) ) ) 

+yshif t ( i ) , ' k ' ) ; 

leftlabel  =  sprintf (' %g ', round (dist (iplot (i) )))  ; 
distlabel  =  sprintf (' %g ', round (dist (iplot (i) ))) ; 
strl  =  sensor ( iplot ( i ),:)  ; 

str2  =  sprintf ( 1 %7 . 2e 1 ,  (xxmax ( iplot ( i ) ) ) ) ; 
if  (xxmax ( iplot  ( i )) )  <  0.0001 

str2=sprintf ( 1  %  4 . lfE-6 ' ,  (lE6*xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 

str2=sprintf ( 1 %4 . 2 fE-3 1 ,  ( IE 3* xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.01 

str2=sprintf ( ' %4 . 2  f E-3 ' ,  ( IE 3* xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf ( 1 %4 . lfE-3 ' , (lE3*xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  10.0 

str2=sprintf ( ' %4 . 2  f ' ,  (xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  100.0 

str2=sprintf (  '  %4 . If ' ,  (xxmax (iplot (i) )  )  )  ; 
elseif  (xxmax ( iplot ( i ))  )  >  1000.0 

str2=sprintf ( 1 %4 . If k 1 , (  xxmax (iplot (i))/1000  )); 
else 

str2=sprintf ( '  %4 . If ' ,  (xxmax (iplot (i) ) )  )  ; 
end  %text  loop 
str3  =  char (strl); 
str4  =  char(str2  ,  str3)  ; 
str5  =  cellstr ( str4 ) ; 

text ( tstartpanell-0 .05-0.02,  [y shift (i) ] , str5,  ' Font Size ' , 8,  .  .  . 

' HorizontalAlignment '  ,  1  right ' ) ; 

text ( tendpanell-0 .02,  [y shift (i) +0 . 3] , distlabel,  1  Fonts ize ' , 8,  . . . 
1 HorizontalAlignment 1 , ' right ' ) ; 

onH  9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9- 
trim.  ooooooooooooooooooooooooooooooooooooooooooooooo 

if  (stype==l ) ; 

plottitle  =  sprintf ('MD  2005  --  Record  %g  Shot  #%s  %s 
shotnum, source2,  srcloc)  ; 

else 

plottitle  =  sprintf ('MD  2005  --  Record  %g  %s  %s 
shotnum, source2,  srcloc)  ; 

end 

subplot  ( 1 , 4 , 2 ) ;  hold  on; 

title  (plottitle,  ' Font Size ' , 18 ,  ' FontWeight '  ,  ' Bold ' ) ; 
set  (gcf ,  ' PaperUnits ' ,  ' inches ' ) ; 

%set (gcf , ' Paper Posit ion ', [0.5,0.5,7,10]); 

% set (gcf , ' PaperOrientation ' , 'portrait ' ) ; 
set (gcf ,  ' Paper Posit ion ',  [0.5,0.5,10,7] )  ; 
set (gcf, ' PaperOrientation ' , ' landscape ' ) ; 

if (iprint==l) ;  print  -dljetplus;  delete (fl);  end 

return  % - 


plotNZ48.m  -  to  plot  the  pistol  or  blasting  cap  data  on  one  panel 
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This  code  plots  48  channels  of  data  in  a  single  panel  on  a  sheet  of  portrait  oriented 
paper.  This  code  does  not  mute  any  data. 

This  code  calls  MD051abel.m  to  get  the  experimental  information  including 
source  type  and  location,  receiver  location,  receiver  sensitivities,  etc. 

%  PlotNZ48.m  -  48  channels 

%  Need  to  set  file  name  and  read  in  data  first 
%  using  readnz  sbr 
%  Calls  sbrs 

%  MD051abel.m  -  supplies  the  details  of  the  expt 
%  including  sensitivity,  coords  of  each  channel  etc. 

%  Don  Albert  Donald.G.Albert@erdc.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 

iprint-1 ;  %  to  print  and  delete  plot 
iprint  =  0;  %skip  printing 

nchan  =  48; 

if (shotnum>=107  &  shotnum  <=117);  nchan  =  36;  end 
if (shotnum>=118  );  nchan  =  24;  end 

srate  =  1/  deltat;  %number  of  samples  per  second 
t  =  ( 1 : npts ) *deltat  +  delayms/1000 ;  % 
tend  =  npts; 

t start 2=1 ; tend2=tend; t2=t ; 
tstartpanell  =  0.0; 
tstartpanell  =  delayms/1000;  % 
tendpanell  =  tstartpanell  +  t (tend) ; 

%  Default  time  window  (entire  record) 

tt2  =  0.0;  tt3  =  t (tend) ; 

tt2  =  tstartpanell;  tt3  =  t (tend) ;  % 

%  Time  windows  for  individual  records 
%if (recnum==235) ;  tt2  =  0.6;  tt3  =  1.4; 

%elseif (recnum==231 ) ;  tt2  =  2.0;  tt3  =  8.0; 

%%elseif (recnum==231 ) ;  tt2  =  1.19;  tt3  =  1.21;  %Lightning  strike? 
%end 

tstart2  =  f loor (tt2/deltat ) ; tend2  =  round (tt3/deltat) ; 
if  (tstart2  <=  1);  tstart2  =  1;  end 
tstartpanell  =  tt2 ; tendpanell  =  tt3; 

%ttstart  =  0.0;ttend  =  2.0;  %time  window  in  seconds  for  plot 

%tend  =  length (x (:, 1 )) ;  %Plot  all  of  data  =  default 
%t  =  delayms/1000  +  ( 1 : tend) /srate ; 

%tstart2=l; tend2=tend; t2=t ; 


xx  =  [  ] ; 

for  i  =  1 : nchan; xmax (i) =max (abs (x (:, i) )) ;  end 

for  i  =  1 : nchan; xmax (i) =max (abs (x (tstart2 : tend2 , i) )) ;  end 
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if(xmax(i)  ==  0);  xmax(i)  =  1;  end 

%Now  set  in  doplotnz.m  shotnum  =  str2num ( fname ( length ( fname) - 
6 : length ( fname) -4 ) ) ; 

%  Get  geometry,  etc  for  this  experiment 

[sensor, sensortype, rtype, sensorx, sensory, sensorz,  .  .  . 

srcloc,  sourcetype, stype, srcx, srcy, srcz, sourcesize,  .  .  . 
engf act , dist , sensor2 , source2 ]  =  MD051abel 
(x, shotnum, nchan) ; 

if (shotnum>=112  &  shotnum  <=1 17);  nchan  =  36;  end 
if (shotnum>=118  );  nchan  =  24;  end 

%if (stype==9) ;  nchan=24;  end  %  Pistol  shots,  B&K  mikes  only 

%  Conversion  factor  from  voltage  to  physical  units  (Pa  or  m/s) 
%xgain  =  descalingf act . / (stackcount . *engf act ) ;  %recorder  units  to 
eng  units 
for  i  =  1: nchan; 

if (stackcount (i) ==0) ;  stackcount (i) =1 ;  end 
%  Convert  recorder  units  to  eng  units 

xgain(i)  =  descalingf act ( i )./( stackcount ( i ). *engf act ( i )) ; 
xx ( : , i )  =  x ( : , i )  *  xgain ( i )  ; 

%  xxmax(i)  =  max (abs (xx) ) ;  %abs  max  for  plotting  -  entire 
trace 

xxmax(i)  =  max (abs (xx (tstart2 : tend2 , i) ) ) ;  %abs  max  for  plot¬ 
ting 
end 

[aa  bb]  =  max (xx) ;  %set  plot  windows  because  manual  TB 
[cc  dd]  =  min (xx) ;  %for  output  file 

tmaxl  =  delayms/1000  +  (bb(l)  ) /srate;  %P-l-0  times  in  seconds 
tmax21  =  delayms/1000  +  (bb (21 )) /srate;  %P-100-0 
tstartl=  floor (tmaxl  *10) /10; 
tstart21=  floor (tmax21*10 ) /10 ; 

%tstartpanell  =  0.0; 

%tendpanell  =  tstartpanell  +  t (tend) ; 

textmax=max (max (t) ) +max (max (t) ) /30; 
textmin=min ( t ) -max (max ( t ) ) / 3  0 ; 

nametext  =  [fNZ  Data  -  1  panel'];  % - FIGURE  -  1  Panel 

fl  =  f igure (' Name ', nametext ) ; 

set (gcf , ' Units ' , ' inches ' , ' PaperUnits ' , ' inches ' ) ; 

%  Pressure  sensors  -  Channels  1-12 
%  iplot=12 : -1 : 1 ;  %Channels  1-12 

%iplot  =  [ 12 : -1 : 1  16  15  14  13  29:36  24 : -1 : 17  25:28  ]; 
iplot  =  [24 : -1 : 1  ] ; 
iplot  =  [13:24  12:-1:1]  ; 

iplot  =  [ 12 : -1 : 1  24:-l:13  36:-l:25  48:-l:37  ]; 

iplot  =  [48 : -1 : 37  36:-l:25  24:-l:13  12: -1:1  ]; 

iplot  =  [37:48  25:36  13:24  1:12  ]; 

%  Switch  channels  close  to  source  for  bare  ground  shots: 
if  ( srcloc ( 1 )  == ' B ') ; 
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iplot  =  [37:48  25:3 
end 

if ( stype==9 ) ;  iplot  = 
mikes  only 

if  (nchan==36) ;  iplot  : 
B&K  mikes  only 
if  (nchan==2 4 ) ;  iplot  : 
mikes  only 


1:12  13:24  ] ; 

[37:48  25:36  ] ;  end  : 

[25:36  13:24  1:12 

[  13:24  1:12  ] ;  end 


Pistol  shots,  B&K 
;  end  %  Pistol  shots, 
%  Pistol  shots,  B&K 


nplot  =  length (iplot) ; 


%  Set  axis  params  here;  depends  on  nplot  -  t2 
yshift  =  (l:nplot)*2  -  2; 

%ax  =  [ 1 2  ( 1 )  max(t2)  -1  2*nplot] ; 

ax  =  [ tstartpanell-0 . 05  tendpanell  -1  2*nplot] ; 

%ax  =  [13-1  2*nplot] ; 
axis  (ax) ; 

axl  =  gca;  set  (axl ,  1 YTick ’ ,  [ ] ) ;  set (axl ,  ' Box ' ,  ' on ' ) ; 
xlabel ( ’ Time,  sec') 


%  Write  trace  data  to  text  file 
fid7  =  fopen ( 1 Plot48Traces .m* ,  1  a  * ) ; 
fprintf ( f id7 , 1 MD  2005  --  Record  %g  %s  %s 
\n 1 , shotnum, source2 , srcloc) ; 
fprintf (fid7, . . . 

1  i,ii  engfact  xgain  xxmax  mV  =  xxmax/engf act 

\n'); 

for  i=nplot:-l:l 
ii  =  iplot  (i) ; 

fprintf  (fid7,  *%2g  %2g  %10.2e  %10.2e  %10.2e  %10.2e  %s  %s  \nf,... 
i,ii,  engfact  (ii),  xgain  (ii),  xxmax (ii), 
xxmax ( ii ) / engfact ( ii ) , . . . 

sensor(ii, :) ,sensor2 (ii, 1 : 16)  ) ; 

%  [iplot;  engfact ( iplot ) ;  xxmax (iplot) ]  * 
end 

fprintf ( fid7 , '  \nf); 
fprintf ( fid7 , '  \nf); 
f close ( f id7 ) ; 


%  plot  data  and  label 

hold  on; 

for  i=l : nplot 

plot  (t,  (xx (1 : length (t) , iplot (i) ) /xxmax (iplot (i) ) ) 
tyshift (i) ,  1 kf )  ; 

leftlabel  =  sprintf (’ %g round (dist (iplot (i) ))) ; 

strl  =  sensor ( iplot ( i ),:)  ; 

str2  =  sprintf ( ' %7 . 2e 1 , (xxmax ( iplot ( i ) ) ) ) ; 

if  (xxmax ( iplot ( i )) )  <  0.0001 

str2=sprintf ( 1  %  4 . lfE-6 f ,  (lE6^xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 

str2=sprintf ( 1 %4 . 2  fE-3 1 ,  ( IE 3* xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.01 

str2=sprintf ( 1 %4 . 2  f E-3 1 ,  ( lE3^xxmax ( iplot ( i ) ) ) )  ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf ( 1  %  4 . lfE-3 f ,  (lE3*xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot  ( i )) )  <  10.0 
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str2=sprintf ( 1 %4 . 2 f ' ,  (xxmax ( iplot ( i ) ) )  )  ; 
elseif  (xxmax ( iplot ( i )) )  <  100.0 

str2=sprintf ( ' %4 . If  1 ,  (xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  >  1000.0 

str2=sprintf ( '  %4 . If  k ' , (  xxmax (iplot (i))/1000  )); 
else 

str2=sprintf ( ' %4 . If ' ,  (xxmax (iplot (i) ) )  )  ; 
end  %text  loop 

%skip  all  scaling,  just  write  peak 
%  str2=sprintf ( ' %4 . le ' , (xxmax (iplot (i) ) ) ) ; 
str3  =  char(strl); 

%str4  =  char(str3  ,  str2);%Two  lines 
str4  =  [str3  str2];  %Concat  on  one  line 

str4  =  [str3  ] ;  %Concat  on  one  line 
str5  =  cellstr (str4) ; 

text ( tendpanell  +  0 .02,  [y shift (i) ] , str5,  ' Font Size  1 , 8,  .  .  . 

' HorizontalAlignment 1 ,  ' left '  )  ; 

%if (find (iplot (i) = [12  11])) 
lef tlabel02  =  [leftlabel  str2]; 

text  ( tstartpanell-0 . 05- 
0.02,  [y shift (i) ] , leftlabel02,  1  Font Size '  ,  8,  .  .  . 

1 HorizontalAlignment 1 ,  ' right 1 ) ; 

end 

%end 

if  (stype==l ) ; 

plottitle  =  sprintf('MD  2005  --  Record  %g  Shot  #%s  %s 
shotnum, source2,  srcloc)  ; 

else 

plottitle  =  sprintf('MD  2005  --  Record  %g  %s  %s 
shotnum, source2, srcloc)  ; 

end 

title  (plottitle,  ' Font Size ' , 18 ,  ' FontWeight 1 ,  ' Bold 1 ) ; 

set (gcf , 1 PaperUnits ' , ' inches ' ) ; 

set (gcf, ' Paper Posit ion ', [0.5,0.5,7,10]); 

set  (gcf ,  ' PaperOrientation 1 ,  'portrait 1 ) ; 

if (iprint==l) ;  print  -dljetplus;  delete (fl);  end 


return  % 


MD051abel.m  -  to  get  sensor  geometry,  sensitivities,  etc 

This  function  is  used  to  provide  information  about  the  experiment  to  the  plot  pro¬ 
grams.  It  is  called  by  plotNZ48FourCol.m  and  by  plotNZ48.m.  This  function 
calls  words.m  to  parse  strings  into  individual  words. 

For  each  record  number,  the  function  returns  the  source  type,  size,  and  location  as 
well  as  an  array  of  sensor  information  including  the  locations,  types,  and  sensitivi¬ 
ties  that  can  be  used  to  convert  sensor  voltages  to  physical  units.  The  function 
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tracks  changes  that  were  made  to  the  sensor  array  during  the  test  (when  sensors 
failed  and  were  replaced). 


function  [sensor, sensortype, rtype, sensorx, sensory, sensorz, . . . 

srcloc, sourcetype, stype, srcx, srcy, srcz, sourcesize, . . . 
engf act , dist , sensor2 , source2 ]  =  MD051abel 
(x, shotnum, nchan) 

%  Fujitsu  C:\aaaNZData\PlotATC05\MD051abel.m 
%  Win2k  D:\FY04Sydet\FY04LargeExplosionTest\MD04\MD041abel.m 
%  Returns  channel  ids  for  NZ  seismograph  used  at  MD  2004  test 
%  calls  words. m  -  breaks  string  into  individual  words 

%  Don  Albert  dalbert@crrel.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 

nchan  =  48; 


%  Sensor  information 
%  Sensor  coordinates 


sensorx  =  ( 1 : nchan) *0 ;  sensory  =  ( 1 : nchan) *0 ;  sensorz  = 
( 1 : nchan) *0 ; 


sensorx  =  [92,92,92,92,  62,62,62,62,  32,32,22,12,  ... 

92,92,92,92,  62,62,62,62,  32,32,22,12,  ... 

240,240,240,240,  180,180,180,180,  120,120,120,120,  ... 

400,400,400,400,  370,370,370,370,  340,340,340,340]; 
sensory  =  [-15,-15,-15,-15,  -15,-15,-15,-15,  -15,-15,-15,-15,  ... 

15,15,15,15,  15,15,15,15,  15,15,15,15,  ... 


0,0, 0,0,  0,0, 0,0,  0,0, 0,0,  0,0, 0,0,  0,0, 0,0,  0, 0,0,0]; 


heightl5  =  [3,8,10,11,  16,20,22,23,  28,32,36,  40,44,48  ]; 
sensorz (heightl5)  =  1.5; 

%  Sensortype 

for  i  =  1  inchan;  sensortype ( i )  =  ' P 1 ;  end 

senv  =  [  2  6  14  18  26  30  34  38  42  46];  sensortype ( senv) =' V 1 ; 
senr  =  [  1  5  13  17  25  29  33  37  41  45];  sensortype ( senr) =  ' R ' ; 
%sent  =  [  12];  sensortype (sent )='T'; 


RP-92 

VP-92 

PP-92-1 

PP-92-0 

! 

. . .  %l-4 

RP-62 

VP-62 

PP-62-0 

PP-62-1 

! 

.  .  .  %5-8 

PP-32-0 

PP-32-1 

PP-22-1 

PP— 10  —  0 

! 

...  %9-12 

RB-92 

VB-92 

PB-92-0 

PB-92-1 

! 

...  %13-1 6 

RB-62 

VB-62 

PB-62-0 

PB-62-1 

! 

...  %17-2  0 

PB-32-0 

PB-32-1 

PB-22-1 

PB-10-0 

! 

...  %2 1-2  4 

R-240 

V-24  0 

P-240-0 

P-240-1 

! 

...  %25-28 

R-180 

V-180 

P-180-0 

P-180-1 

! 

...  %2  9-32 

R-12  0 

V-120 

P-120-0 

P-120-1 

! 

...  %33-3  6 

R-4  00 

V-4  00 

P-400-0 

P-400-1 

! 

...  %37-4  0 

R-37  0 

V-37  0 

P-370-0 

P-370-1 

! 

...  %4 1-4  4 

R-340 

V-34  0 

P-340-0 

P-340-1 

! 

];  %45- 

sensor  =  words (sss)  ; 

%  Initial  sensor  array,  changes  during  test  listed  below 
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%  Corrected  location  of  15974  on  the  ground,  15972  at  1.5  m 

height 

sss2  = [  ... 

' MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  PCB-10 6B50-6522  PCB- 
10  6B50-6693  '  ...%l-4 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  PCB-102A07-15973  PCB- 
102A07-15971  '  ...%5-8 

' PCB-102A0  9-15  968  PCB-102A0 9-15 969  PCB-102A0  6-15  964  PCB- 

102A09-15962  1  ...%9-12 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  PCB-106B50-3259  PCB- 
10  6B50-6695  '  ...%13-16 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  PCB-102A07-15974  PCB- 
102A09-12102  1  . . . %17-20 

' PCB-102A0 9-131 80  PCB-102A0 9-15 970  PCB-102A0 6-15 966  PCB- 

102A06-15965  '  ...%21-24 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  B&K-4938-6  B&K-4938-5 

1  . . . %25-28 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  B&K-4938-4  B&K-4938-3 

'  .  . . %2  9-32 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  B&K-4138-2  B&K-4938-1 

1  .  .  . %33-3  6 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  B&K-4165-6  B&K-4965-5 

'  . . . %37-40 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  B&K-4165-10  B&K-4165-9 

1  .  .  . %4 1-4  4 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  B&K-4165-8  B&K-4165-7 

1  ];  %45-48 

sensor2  =  words (sss2); 

%  rtype  1  =  Microphone,  2  =  PCB,  5=  Vert,  6  =  Horiz,  7  =  Accel 
rtype  =[6522  6522  2222  6522  6522  2222 

6511  6511  6511  6511  6511  6511]; 

%  engfact  =  units  per  mV 

%  mV  =  data_value*descal_f act/stack_count 

%  phys_units  =  data_value*descale_f act  /  ( stack_count*engf act ) 

%Starting  values : 
engfact ( 1 : nchan)  =  1.0; 

engfact ([27  28  31  32  35  36  39  40  43  44  47  48])  =  -1.0;  %B&K  po¬ 
larity 

%  Calibrated  values:  mV  per  m/s  and  mV/Pa 

%  Calibrated  values:  --  B&K  field  cal  and  PCB  cals  at  CRREL 


(mV/Pa) 
engfact  = 
32.2*1000 

-32.2*1000 

0.0668 

0.0636 

V 

1 

1 — 1 

o\° 

32.2*1000 

-32.2*1000 

11.5  6e-3 

12 . 1 6e-3 

.  .  .  %5-8 

5 . 77e-3 

5. 60e-3 

1 . 23e-3 

1 . 21e-3 

...  %9-12 

32.2*1000 

-32.2*1000 

0.0778 

0.0630 

...  %13-1 6 

32.2*1000 

-32.2*1000 

11 . 91e-3 

5 . 14e-3 

...  %17-20 

5 . 42e-3 

5 . 8  6e-3 

1 . 22e-3 

1 .22e-3 

...  %2 1-2  4 

32.2*1000 

-32.2*1000 

-1.13/10 

-1.24/10 

...  %25-28 

32.2*1000 

-32.2*1000 

-1.22/10 

-1.11/10 

...  %2  9-32 

32.2*1000 

-32.2*1000 

-1.12/10 

-1.05/10 

...  %33-3  6 
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32 .2*1000 

-32 .2*1000 

-49.1/100 

-47.5/100  . 

.  .  %37-4 0 

32 .2*1000 

-32 .2*1000 

-48 . 5/100 

-49.6/100  . 

.  .  %4 1-4  4 

32 .2*1000 

-32 .2*1000 

-61.2/100 

-57 . 4/100 

];  %45-48 

%  Changes  to  sensor  array  during 
if ( shotnum>=4 0 ) ;  sensor2  (2 1 ,  : )  = 


engf act (21 ) 
if  (shotnum>=43 ) 
engf act (21 ) 
if  (shotnum>=47 ) 
engf act ( 12 ) 
if  (shotnum>=47 ) 
engf act ( 11 ) 


tests 

' PCB-102A0 9-12 101 


=  4.88e-3;  end 
sensor2 (21 , : ) 

=  1.16e-3;  end 
sensor2  ( 12 ,  : ) 

=  1 . 2  8e-3 ;  end 
sensor2  (11,  : ) 

=  1.26e-3;  end 


=  ' PCB-102A0  6-1 67  84 


=  ' PCB-102A0  6-15503 


' PCB-102A0  6-15961 


%  Source  information  - 

%  Default  source  info 
srcloc  =  1 SP1 1 ; 

srclocnum  =  1;  srcx  =  0.0;  srcy  =  0.0;  srcz  =  1.5; 

%  stype  1  =  C4,  6  =  50  cal,  9  =  45  cal,  10  =  noise 
%  stype  3  =  Claymore,  2  =  DAM 

%  stype  12  =  Hammer,  13  =  Cap,  14  =  Calibration 
sourcetype  =  'C4';  stype  =  1;  sourcesize  =  0; 
source2  =  1  Jungle ' ; 

%  MD05:  All  shots  were  1  brick  of  C4 

sourcetype  =  '  C4';  stype  =  1;  sourcesize  =  1; 


%  C4  shots 

if (  shotnum 
shotnum 
srcloc  = 

end 

if  (  shotnum 
shotnum 
srcloc  ~ 

end 

if  (  shotnum 
shotnum 
srcloc  = 

end 

if  (  shotnum 
shotnum 
srcloc  = 

end 

if  (  shotnum 
shotnum 
srcloc  = 

end 

if  (  shotnum 
shotnum 
srcloc  = 

end 

if  (  shotnum 
srcloc  = 

end 

if  (  shotnum 
srcloc  = 

end 


==38  |  shotnum  ==40  |  shotnum  ==75  |  ... 

==77  |  shotnum  ==79) 

' B+2+1.5';  srcx  =  2.0;  srcy  =  15.0;  srcz  =  1.5; 


==39  |  shotnum  ==42  |  shotnum  ==76  |  ... 

==78  |  shotnum  ==80) 

1  P+2  +  1.5 f;  srcx  =  2.0;  srcy  =-15.0;  srcz  =  1.5; 


==37  |  shotnum  ==43 |  shotnum  ==45  |  ... 

==81  |  shotnum  ==83  |  shotnum  ==85) 

' B-5+3';  srcx  =  -5.0;  srcy  =  15.0;  srcz  =  3.0; 

==44  |  shotnum  ==46  |  shotnum  ==82  |  ... 

==84  |  shotnum  ==86) 

'P-5+3';  srcx  =  -5.0;  srcy  =-15.0;  srcz  =  3.0; 


==69  |  shotnum  ==71  |  ... 

==87  | shotnum  ==89  |  shotnum  ==91) 

'B-5+1 ' ;  srcx  =  -5.0;  srcy  =  15.0;  srcz  =  1.0; 


==70  |  shotnum  ==72  |  ... 

==88  | shotnum  ==90  |  shotnum  ==92) 

'P-5+1';  srcx  =  -5.0;  srcy  =-15.0;  srcz  =  1.0; 


==73  ) 

' B-5+1. 5';  srcx  =  -5.0;  srcy  =  15.0;  srcz  =  1.5; 


==74  ) 

'P-5+1. 5';  srcx  =  -5.0;  srcy  =-15.0;  srcz  =  1.5; 
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if ( (shotnum  >=37  &  shotnum  <=  40));  source2=num2str ( shotnum-3 6 ) ; 
end 

if ( (shotnum  >=42  &  shotnum  <=  46));  source2=num2str (shotnum-37 ) ; 
end 

if ( (shotnum  >=69  &  shotnum  <=  92));  source2=num2str (shotnum-59)  ; 
end 

%  Other  source  info 

if (shotnum  <=8  |  (shotnum  >=12  &  shotnum  <=  35)  |  ... 

(shotnum  >=  93  &  shotnum  <=  104)  |  ... 

(shotnum  >=112  &  shotnum  <=  147)  ) 

sourcetype  =  ' Calibration '; stype  =  14;  srcz  =  0.0; 
source2= 1  Cal 1 ; 
end 

if (shotnum  ==106  |  shotnum  ==107  |  shotnum  ==108) 

sourcetype  =  'Blasting  Cap'; stype  =  15;  source2= ' Cap ' ; 
srcloc  =  'P+2+2';  srcx  =  2.0;  srcy  =-15.0;  srcz  =  2.0; 

end 

if (shotnum  ==109  |  shotnum  ==110  |  shotnum  ==111) 

sourcetype  =  'Blasting  Cap'; stype  =  15;  source2= ' Cap ' ; 
srcloc  =  'BP+2+2';  srcx  =  2.0;  srcy  =  15.0;  srcz  =  2.0; 

end 

if (  shotnum  ==36  |  shotnum  ==41  |  shotnum  ==105  |  shotnum  ==47  | 

shotnum  ==68  |  shotnum  ==999  |  shotnum  ==999  |  shotnum 
==999)  ; 

sourcetype  =  ' Noise '; stype  =  10;  srcz  =  0.0;  source2= ' Noise ' ; 
srcloc  =  '  ' ; 

end 

if (shotnum  >=8  &  shotnum  <=11) 

sourcetype  =  ' Vehicle '; stype  =  15;  source2= ' Vehicle ' ; 
srcloc  =  'Road';  srcx  =  240.0;  srcy  =  100.0;  srcz  =  3.0; 

end 

%  Blank  pistol  shots  along  the  sensor  array 
if (  (shotnum  >=48  &  shotnum  <= 6 7 ) ) 

sourcetype  =  '45  cal  blank';  stype  =  9;  source2= ' Pistol ' ; 

sourcesize  =  0;  srcz  =  1.0; 

if (  (shotnum  >=48  &  shotnum  <=52) ) 

srcx  =  280.0;  srcloc  =  '280  m';  srclocnum  =  7; 

end 

if (  (shotnum  >=53  &  shotnum  <=5 7 ) ) 

srcx  =  310.0;  srcloc  =  '310  m';  srclocnum  =  8; 

end 

if (  (shotnum  >=58  &  shotnum  <= 6 2 ) ) 

srcx  =  430.0;  srcloc  =  '430  m';  srclocnum  =  9; 

end 

if (  (shotnum  >=63  &  shotnum  <= 6 7 ) ) 

srcx  =  280.0;  srcloc  =  '280  m';  srclocnum  =  7; 

end 

end 

%  Source  to  sensor  distances  - 

srcx2  =  srcx*srcx;  srcy2  =  srcy*srcy;  srcz2  =  srcz*srcz; 
for  i=l : nchan 

dist(i)  =  ( (sensorx (i) -srcx) . A2  +  ... 
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( sensory ( i ) -srcy) . A2  +  ... 

( sensor z ( i ) -srcz ). A2  )  .A0.5; 

end 

return  % - 
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Appendix  B.  StrataVisor  NZ  Recorder  Specifications 


Strata  Vis  or™ 


□  Professional,  ruggedized  24  bit  seismic 
recorder  suitable  for  reflection,  refraction, 
downhole,  VSP  or  small  marine  surveys. 
Intelligent  self-trigger  enables  vibration  and 
earthquake  monitoring. 

□  Houses  from  3  to  64  channels.  Expands 
seamlessly  up  to  1000  channets  by 
connecting  Geode  in-field  distributed 
modules, 

□  Engineered  for  rugged  field  use:  reliable  in 
harsh  environments,  survives  shock,  humidity 
and  dust.  Passes  MIL  810E  vibration  spec. 

□  20  kHz  bandwidth  (0.02  jjs  to  16  ms  sampling) 
for  ultra-high  resolution  engineering  surveys 
or  low  frequencies  for  earthquake  monitoring. 

□  Brilliant  fuEE-sun-visihle  color  screen  and  built- 
in  plotter  -  available  also  as  rugged  field 
computer  without  seismic  channels. 

□  Built-in  geophone  and  line  testing,  full 
waveform  noise  monitor.  Optional  automated 
internal  in-field  instrument  testing  or 
detachable  lab  quality  externaF  test  oscillator. 

□  High-speed  hardware  correlator  means  no 
delays  using  Vibroseis  or  pseudo-random 
Mini-Sosie  sources. 

□  Powerful  reflection,  refraction  and  cross  hole 
tomography  software: 

■  M  od  eT  local  geology  a  nd  d  isp  la  y  tra v  el  -ti  m  e  cu  rye  s 
before  surveying 

■  Pick  breaks  on  site  and  view  travel-time  curves  for 
optimum  shot  positioning 

■  Display  in-field  preliminary  cross  section  and  ray 
coverage  plots  to  see  what  you  have  missed 

■  Undertake  a  comprehensive  analysis  back  at  the 
office  and  easily  compare  results  with  alternate 
interpretation  methods 


NZ  Exploration  Seismograph  and  Field  PC 


The  StrataVisor1 NZ  is  a  high-performance  exploration 
seismic  system  in  a  compact,  weatherproof  chassis.  The 
NZ  can  operate  as  a  field  PC,  as  a  stand-alone  seismic 
recorder  with  3  to  64  internal  channels  -  or  expands  easily 
to  larger  channel  systems  by  connecting  Geode  remote 
distributed  modules.  This  flexibility  lets  you  collect  data  for 
all  applications  in  all  environments  -  you  can  even  rent 
extra  channels  when  needed. 


Examine  your  data  at  all  phases  of  acquisition  to  ensure 
data  quality.  Customizable  windows  show  real-time  noise 
monitor,  amplitude  spectra  and  seismic  traces  so  you  see 
problems  instantly.  A  log  file  keeps  track  of  all  parameter 
changes  and  customizable  alarms  alert  you  to  critical 
issues.  You  can  even  do  preliminary  processing  in  the  field 
with  industry-leading  reflection,  refraction  and  tomography- 
software  included  with  every  system. 

The  StrataVisor  NZ  console  includes  a  brilliant  daylight- 
visible  color  screen,  waterproof  keypad  and  built-in  printer. 
Low- power  circuitry  and  a  standby  mode  extend  battery  life 
and  reduce  weight.  A  tape  drive  controller,  Vibroseis  sweep 
generator  and  intelligent  integration  of  GPS  and  other 
survey  data  can  be  included. 

The  StrataVisor  NZ  is  backed  by  a  I  year  parts  and  labor 
warranty.  All  this  from  a  company  with  factory  trained 
service  centers  world  wide  and  over  30  years  of  superior 
support  to  geoscience  professionals. 


■  *  *  * 


*  •  *  * 


hi* 
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Specifications: 


Strata  Visor™  NZ  Exploration  Seismic  Recorder 


Configurations: 

■  Lightweight  field-iugged  PC  with,  no  seismic  rlim-jl-.  for  use  as  field 
computer  or  nigged  in-field  controller  far  Geode  distributed  modules 

■  Integrated  seismic  recorder  -  add  3.  6,  B.  12  or  16  to  64  built-in  channels 
(in  8  channel  increments  : 

■  Connect  multiple  NZ :  s  and  operate  Horn  a  single  keypad 

Runs  Windows™  operating  system  and  includes  all  software  far  controlling 
internal  channeh  and  multiple  lares  of  Geode  modules.  Total  number  of  channels 
limited  only  by  practical  survey  requ  ire  rue  nts . 


Continuous  Recording:  Available  for  vibration  monitoring. 

Delay:  0  to  ICO  sec  in  1  sample  steps. 

Auxiliary  Channels:  All  channels  can  be  programmed  as  either  AUX  or  DATA. 

CDP  Roll  Along:  Software  selectable  channels  can  be  robed  through  total 
channels. 


AD  Conversion:  24-  bit  result  using  Ciystal  Semiconductor  sigma-delta 
converters  and  Geometric  s  proprietary  over- sampling . 


Dynamic  Range:  1+4  dE  (system).  110  dB  (instantaneous, 
measured)  at  2  ms.  24  dB. 

Distortion;  0.0005%  @  2  ms.  1.75  to  208  Hz. 

Bandnidth:  1.75  Hr  to  20  kHz.  Low  comer  frequency  option 
available. 

Common  Mode  Rejection:  >  lQOdB  at  -  100  Hz.  36  dB. 

Cross  Talk:  -125  dB  at  23.5  Ha.  24  dB.  2  ms. 

Noise  Floor  0.20  (.lY,  RH  at  2  ms.  36  dB.  1.75  to  20®  Hz. 

Stacking  Trigger  .Accuracy:  1/32  of  sample  interval. 

Maximum  Inpnt  Signal:  2.SY  FF.  0  dB.,  177  mYPP.  24  dB. 

Inp  nt  Imp  ed  ance :  20  kQhm  0.02  |+f. 

Preamplifier  Gains:  Standard  factray  configuration  is  2+  and 
36  db.  selectable  in  software.  Optionally,  can  be  jumpered  for 
software  selectable  12  and  24  dB  or  can  be  jumpered  in  four 
channel  blocks  a.s  a  single  faced  gain  of  0  dB  for  high-voltage  devices. 


Line  Testing:  null  waveform  noise  monitor  displays  real-time  output  from 
geophones.  Tests  geophones  and  locates  shorted  or  broken  cables. 


Optional  Built-In 

Test  Functions 
Digital 

*  CPU  diagnostics 

*  Internal  network  test 

*  Digital  nocuous 

*  Battery  Warning 
Instrument: 

*  Noise 

*  DC  Offset 

*  Gaii  Accuracy 

*  Gain  and  Phase 
Similarity 

*  Distortion. 

*  Cnossfeed 

*  Bandwidth 

+  Timing  Accuracy 
Line: 

*  Geophone  similarity 
+  Noise 

*  Disconnected  phones 


Instrument  Tests:  Available  either  internally  or  is  detachable 
external  system.  Test  scripts  are  customizable  to  meet  your  clients 
needs. 

Data  Formats:  SEG-2.  SEG-D  and  SEG-Y. 

System  Software:  Runs  under  Windows™  operating  system. 

Uses  Geometries  MGOS  software  to  control  acquisition  on  internal 
channels  and  up  to  4  hues  of  external  channels  housed  in  Geode 
distributed  modules.  System  includes  the  following  2lkl  paity 
applicahons.  software: 

■  S3FQC  delay  time  refraction  software  from  Rimnock 

Geophysics 

■  SeisImaeeT  left-action  modelinE  and  analysis  software  horn 

OYG 

■  WmSeis-Ti.ubo  reflection  software  fiom  the  Kansas  Geological 

Survey 

Redaction  software  packages  mstalledon  mstnunent  aie 
coufigined  for  m-fleld  analysis  and  may  require  attaching  a 
keyboard  and  mouse.  Full  desktop  versions  of  these  packages  are 
available  either  faom  Geometries. 

Data  Storage:  Stores  data  on  internal  hard  drive  or  controls 


multiple  SCSI  tape  devices.  Records  in  SEG2.  SEGY  or  SEGD. 


.Anti-dhaiS  Filters:  -3  dB  aft  B3° :  of  Nyquist  down  93  dE. 

.Acquisition  and  Display  Filters: 

Low  Cut:  OUT.  10,  15, 25,  35:  50,  70:  100, 140,  20O:  2E0:.  400  Hz.  24  or  48 
dB' octave.  Butterworth. 

Notch:  50,  60, 150,  180  Hz  and  OUT.  with  the  50  dB  rejection  bandwidth 
2%  of  center  frequency. 

High  Cut:  OUT,  250, 500  or  1000  Hz:  24  or  4®  dB.  octave. 

Customer  filter  fiequencies  available  as  an  option. 

Sample  Interval:  0.02, 0.03125, 0.0625, 0.125, 0.25, 0.5, 1.0, 2.0, 4.0,  S  O,  16  0 
ms. 

Correlation  :  Built-in  high-speed  hardware  correlator  for  Yibnoseis.  Optional 
pile:  conditioning  for  acquisition  of  pseudo-random  (MiniSosie)  somces. 

Record  Length:  16.384  samples  standard.  65.536  samples  optional. 

Fre- trigger  Data:  Up  to  full  record  length. 

Intelligent  Self-Trigger:  Available  for  earthquake  and  vibration  monitoring. 


Plotters:  Built-in  4"  thermal  plotter.  Drives  a  variety  of  ITT  compatible  plotters 
including  Pitntrex  4.  3  and  12  inch  continuous  thermal  plotters. 

Triggering:  Positive,  negative  or  contact  closure,  software  adjustable  threshold. 

Power:  30Wplus  3. 65  W  dbannel  during  acquisition.  Standby  mode  reduces 
channel  power  consumption  by  70%.  Requires  external  12V  sipply. 

Emironmenta]:  Boots,  fiom  -5°C  to  4^C.  Operates  from— 5*0  to  4D"C. 

Extended  temperature  version  available  to  -6G°C .  Operates  in  a  light  rain, 
watertight  with  cover  closed.  Passes  MILS  1  DEB  vibration  test 

Physical: 

Field  PC  with  no  seismic  channels:  13.: "Lx  18"Wx  13“D  (27crcLx46cm 
w  x  46  cm  D).  weighs  27  Lb  ( 12.3  kg; 

Seismic  recorder  with  3-64  internal  channels :  10.5  "L  x  1S"W  x  21  ‘:D 
(26.7cm Lx  34  cm  W  x  33  cm  D),  weighs  38  lb  (13  !kg) 

Wairaatj-:  One  year- parts  and  labor-.  Extended  warranty  available. 
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